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Chapter 1 



Introduction 



Quantum world opens up several puzzling aspects that are not amenable to classical 
intuition. Correlation exhibited by subsystems of a composite quantum state is one such 
striking feature and has been a source of philosophical debates - following the famous 
Einstein-Podolsky-Rosen discussion on the foundational aspects of quantum theory. It is 
now well established that entangled states play a crucial role in the modern quantum infor- 
mation science, including quantum cryptography QJ, quantum communication and quan- 
tum computation El [3] |4j] . 

Considerable interest has been evinced recently j5lloll7l[^l^[T0l[TTl[T2l[T3l[T4l[T5l 
in producing, controlling and manipulating entangled multiqubit systems due to the possi- 
bility of applications in atomic interferometry lTT6l fTTIl . high precession atomic clocks ||T8l . 
quantum computation and quantum information processing J2j. Multiqubit systems, which 
are symmetric under permutation of the particles, allow for an elegant description in terms 
of collective variables of the system. Specifically, if we have N qubits, each qubit may be 
represented as a spin-^ system and theoretical analysis in terms of collective spin operator 

N 

J = \ ® a (® a denote the Pauli spin operator of the a th qubit), leads to reduction 

a=l 

of the dimension of the Hilbert space from 2^ to (N + 1), when the multiqubit system 
respects exchange symmetry. A large number of experimentally relevant multiqubit states 
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exhibit symmetry under interchange of qubits, facilitating a significant simplification in 
understanding the properties of the physical system. While complete characterization of 
multiqubit entanglement still remains a major task, collective behavior such as spin squeez- 
ing ISIEUIMEMMMMMMMMMMMMMMMM, exhibited by 

multiqubit systems, has been proposed as a signature of quantum correlation between the 
atoms. A connection between spin squeezing and the nature of quantum entanglement has 
been explored ||2"9l |30l and it is shown that the presence of spin squeezing essentially re- 
flects pairwise entanglement. However, it is important to realize that spin squeezing serves 
only as a sufficient condition - not a necessary one - for pairwise entanglement. There 
will still be pairwise correlated states, which do not exhibit spin squeezing. In a class of 
symmetric multiqubit states it has been shown ll30l that spin-squeezing and pairwise 
entanglement imply each other. Questions like"Are there any other collective signatures 
of pairwise entanglement? " are still being investigated. Recently, inequalities generaliz- 
ing the concept of spin squeezing have been derived [31]. These inequalities are shown 
to provide necessary and sufficient conditions for pairwise entanglement and three-party 
entanglement in symmetric iV-qubit states. 

In this thesis, we have addressed the problem of characterizing pairwise entangle- 
ment in symmetric multiqubit systems in terms of two qubit local invariants ll32l l33l l34l . 
This is important because quantum entanglement reflects itself through non-local correla- 
tions among the subsystems of a quantum system. These non-local properties remain unal- 
tered by local manipulations on the subsystems and provide a characterization of quantum 
entanglement. 

Two composite quantum states p\ and p2 are said to be equally entangled if they are 
related to each other through local unitary operations, which merely imply a choice of 
bases in the spaces of the subsystems. One may define a polynomial invariant, which is by 
definition any real valued function of density operators, taking the same value for equally 
entangled density operators p. 
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Basic issues of importance would then be 

• to find complete set of polynomial entanglement invariants which assume identical 
values for density operators related to each other through local unitary operators. 

• decide whether the set of separable states can be described in terms of a polynomial 
invariant /, such that f(p) > is equivalent to separability ll35l . 

In this context, Y. Makhlin |f36| has studied the entanglement invariants of an arbitrary 
mixed state of two-qubits and has identified a complete set of 18 local invariants character- 
izing the system. A set of 8 polynomial invariants has been identified in the case of pure 
three qubit states |[37| . Linden et. al. |38| have outlined a general prescription to identify the 
invariants associated with a multi particle system Q. Here, we focus on constructing a com- 
plete set of local invariants characterizing symmetric two qubit systems and analyzing the 
pairwise entanglement properties like collective spin squeezing - exhibited by multiqubits - 
in terms of two qubit entanglement invariants. 

A brief Chapter wise summary of the thesis is given below. 

Chapter 2: Symmetric two qubit local invariants 

For an arbitrary two qubit mixed state, Makhlin [36] has proposed a complete set of 
18 polynomial invariants. In this Chapter, we show that the number of invariants reduces 
from 18 to 6 in the case of symmetric two qubit states owing to the exchange symmetry. 
We quantify entanglement in symmetric two qubit states in terms of these complete set of 
six invariants. More specifically, we prove that the negative values of some of the invariants 
serve as signatures of quantum entanglement in symmetric two qubit states. This leads us to 
identify sufficient conditions for non-separability in terms of entanglement invariants ||32l . 
Further, these conditions on invariants are shown here to be both necessary and sufficient 
for entanglement in a class of symmetric two qubit states. 

'However, separability properties of two qubit states in terms of the local invariants, is not investigated in 
Ref. [36|[32][3U 
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Chapter 3: Characterization of pairwise entanglement in symmetric multiqubit sys- 
tems 

As discussed in Chapter 2, some of the symmetric two qubit invariants reflect nonsep- 
arability ll32l . In this Chapter, we focus on the characterization and classification of pair- 
wise entanglement in symmetric multi-qubit systems, via local invariants associated with 
a random pair of qubits drawn from the collective systems. In other words, we investigate 
collective signatures of pairwise entanglement in symmetric N-qubit states as implied by 
the associated non-positive values of the two qubit invariants. More specifically, we iden- 
tify here that a symmetric multi-qubit system is spin squeezed iff one of the entanglement 
invariant is negative. An explicit classification, based on the structure of local invariants for 
pairwise entanglement in symmetric A r -qubit states is given ll33l . We show that our char- 
acterization gets related to the generalized spin squeezing inequalities of Korbicz et. al OTTl . 

Chapter 4: Analysis of few dynamical models 

In the light of our characterization of pairwise entanglement in symmetric multiqubit 
states discussed in Chapter 3, we analyze some of the experimentally relevant N qubit 
permutation symmetric states and explicitly demonstrate the non-separability of such states 
as exhibited through two qubit local invariants. In particular, we evaluate the two qubit 
local invariants and hence discuss the collective pairwise entanglement properties in the 
following multiqubit states: 

1. Dicke states 0911401 

2. Kitagawa-Ueda state generated by one axis twisting Hamiltonian |fT9l 

3. Atomic squeezed states |[20l . 
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Chapter 5: Necessary and sufficient criterion in symmetric two qubit states 

Continuous variable systems (CV) [41] i.e., systems associated with infinite dimen- 
sional spaces are a focus of interest and attention due to their practical relevance in appli- 
cations to quantum optics and quantum information science. Moreover two mode Gaussian 
states, a special class of CV systems provide a clean framework for the investigation of 
nonlocal correlations. Consequently, most of the results on CV entanglement have been 
obtained for Gaussian states. Entanglement for two-mode Gaussian states is completely 
captured in its covariance matrix . It is desirable to look for an analogous covariance matrix 
pattern in finite dimensional systems - in particular in multiqubits. 

In this Chapter, we identify such a structural parallelism ll34l between continuous 
variable states [423 and symmetric two qubit systems by constructing covariance matrix of 
the latter. Pairwise entanglement between any two qubits of a symmetric N qubit state is 
shown to be completely characterized by the off-diagonal block of the two qubit covariance 
matrix. We establish the inseparability constraints satisfied by the covariance matrix P4l 
and identify that these are equivalent to the generalized spin squeezing inequalities |f3TTl 
for pairwise entanglement. The interplay between two basic principles viz, the uncertainty 
principle and the nonseparability gets highlighted through the restriction on the covariance 
matrix of a quantum correlated two qubit symmetric state. So, the collective pairwise en- 
tanglement properties of symmetric multiqubit states depends entirely on the off diagonal 
block of the covariance matrix. We further establish an equivalence between the Peres- 
Horodecki ll43l l44l criterion and the negativity of the covariance matrix C showing that our 
condition is both necessary and sufficient for entanglement in symmetric two qubit states. 
We continue to identify the constraints satisfied by the collective correlation matrix V W 
of pairwise entangled symmetric N qubit states. 
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In other words, the local invariant separability condition necessarily implies that 

The symmetric N qubit system is pairwise entangled iff the least eigen value of the 
real symmetric matrix V^ N > + -k SS T is less than N / 4. 

(y( N ) denotes the collective covariance matrix and S corresponds to the collective average 
spin of the symmetric N-qubit system.) 
Chapter 6: Summary 

In this Chapter, we briefly summarize the important results obtained in this thesis. 



Chapter 2 



Symmetric two qubit local invariants 

2.1 Introduction 

Initiated by the celebrated Einstien-Podolsky-Rosen criticism B3I . counterintuitive features 
of quantum correlations have retained the focus for more than seven decades now, and 
quantum entanglement has emerged as an essential ingredient in the rapidly developing area 
of quantum computation and quantum information processing (2l[3j|4l. Characterization 
and quantification of entanglement has been one of the central tasks of quantum information 
theory. In simple terms, a bipartite quantum system is entangled, if it is not separable i.e., 
if the density matrix cannot be expressed as a convex mixture of product states, 

P = y^^Pw Pw ] ® Pw ] where < p w < 1 and = 1. (2.1) 

w w 

Here, {pw^} and {pffl} denote a set of density operators associated with quantum systems 
1 and 2. It is a non-trivial task to check whether a given state is expressible as a mixture of 
product states (see Eq. (12.11 )) or not. 

Peres f43 ] has identified that the partial transpose of a separable bipartite state p is pos- 
itive definite (See Appendix |B]) and therefore negative eigenvalues of a partially transposed 
density matrix imply non-separability of a quantum state. Further, Horodecki et. al H4l 

7 
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proved that negativity under partial transpose provides a necessary and sufficient condition 
for quantum entanglement in 2 g) 2 and 2 (g) 3 systems only. 

It is possible to quantify the amount of entanglement in a bipartite pure state \tp) 
through the von Neumann entropy of either of the two subsystems ll46l 



where, 



and 



E{$) = -Tr(p A \og 2 p A ) = -Tr(p B \og 2 p B ) (2.2) 



p A = TT B (\i>)(i>\) 



p B = Tr A (\i>)(i>\) 

denote subsystem density matrices. E(ip) is referred to as Entropy of entanglement. 
linin^oo nE{ijj) gives the number of maximally entangled states that can be formed with n 
copies of IV 7 ) 5 in the asymptotic limit. 

The entanglement of formation of a mixed bipartite state p is defined as the minimum 
average entanglement 

E(p) = mm (2.3) 

i 

where corresponds to all possible decompositions of the state through 

i 

Entanglement of formation E{p) reduces to entropy of entanglement E(ip) in the case of 
pure states and is zero iff the state is separable. 

An explicit analytical expression for the entanglement of formation has been derived 
for an arbitrary pair of qubits B71I481 and is given by: 



BW _Ji±4±^, (2.4) 



Symmetric two qubit local invariants 
where 

h(x) = —x log 2 x — (1 — x) log 2 (l — x) 
In Eq. (12.41 ). C, the Concurrence ll48l is given by 

C = max(0, Ai — A2 — A3 — A4), 

with Af , A|, A§, A| denoting the eigenvalues of 

p(o 2 <8> o- 2 )p*(cr 2 <£> «t 2 ) 

in the decreasing order. Here, 



0"2 



V 



-i 

1 



(2.5) 



is the standard Pauli matrix and p is expressed in the standard two qubit basis set 

{|0i 2 ),|0i 1 2 ), |li 2 >, |li 1 2 >}- 

The Concurrence C varies from zero to one and is monotonically related to entanglement 
of formation E(p), thus gaining the status of a measure of entanglement on its own 11481 . 

Entanglement properties of a quantum system remain unaltered when the subsystems 
are locally manipulated and two quantum states p\ and p 2 are equally entangled if they are 
related to each other through local unitary transformations. Non-separability of a quantum 
state may thus be represented through a complete set of local invariants which contains 
functions of the quantum state that remain unchanged by local unitary operations on the 
subsystems. In this Chapter, we investigate a complete set of local invariants for arbitrary 
symmetric two qubit states. 

We identify that a set of six invariants is sufficient to characterize a symmetric two 
qubit system, provided the average spin | (a) \ of the qubits is non-zero. If | (a) \ = 0, only 
two entanglement invariants represent the nonseparability of the system. 



Symmetric two qubit local invariants 10 

We further show that all the invariants associated with separable symmetric systems 
are positive. This allows us to identify criteria for non-separability in terms of invariants. 



2.2 Arbitrary two qubit density matrix 

Density matrix (See Appendix |A]) of an arbitrary two-qubit state in the Hilbert-Schmidt 
space H = C 2 <8) C 2 is given by 



1 / 3 3 3 \ 

- / ® I + Si au + ^2 a 2i n + ^2 tij a u a 2 j , 

\ i=l i=l *,J=1 / 



(2.6) 



where I denotes the 2x2 unit matrix. Here, 



1 

0-1=1 1,0-2 

1 



-i 

1 



, 0- 3 



1 

-1 



(2.7) 



are the standard Pauli spin matrices and 



<y\i = ai <g> I 

(?2i = I® CTi- 



(2.8) 



The average spins of the qubits s = (s\, s 2 , S3) and r = (ri, r 2 , r^) are. given by 



Si = Tr (p an 
n = Tr (p a 2i 



(2.9) 



and the two-qubit correlations are given by, 



Uj = Tr [p(aua 2 j)]. 



(2.10) 
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It is convenient to express the two qubit correlations Uj(i, j = 1, 2, 3), in the form 
of a 3 x 3 matrix as follows: 



ft t t ^ 

111 £l2 tl3 

hi *22 *23 

V *31 *32 *33 J 



(2.11) 



The correlation matrix T is real and in general, nonsymmetric. 

Transformation of state parameters under local unitary operations: 

The 15 parameters {sj, n, = 1,2,3, characterizing two qubit density matrix of 

Eq. (12.61 ) exhibit the following transformation properties under local unitary operations: 



i=i 

3=1 



(2.12) 



4. = Y.O^oftM or T' = OWTO( 2 ) T , (2.13) 
fc,z=i 

where O^, E 50(3,i2) are the 3x3 rotation matrices, uniquely corresponding to 
the 2x2 unitary matrices Ui G SU(2). The above transformation properties, facilitate 
the construction of polynomial functions of state parameters {sj, r^, tij} which remain 
invariant ll36l under local operations on individual qubits. We devote the next section for 
discussion of a complete set of local invariants associated with an arbitrary two qubit density 
matrix- which was proposed by Makhlin P6l . 
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2.3 Local invariants of an arbitrary two qubit system 



As has been emphasized earlier, genuine nonlocal properties should be described in terms 
of physical quantities that are invariant under local unitary operations. Makhlin P6l in- 
vestigated such local invariant properties of mixed states of two-qubit system. Two density 
matrices p\ and p 2 are called locally equivalent if one can be transformed into the other by 
local operations 

P2 = {Ui ® U 2 )pi(U 1 <g) u 2 ) ] . 

A useful tool for verification of local equivalence of two states is a complete set of invariants 
that distinguishes all inequivalent states: 

If each invariant from the set has equal values on two states pi, p 2 , their local equiv- 
alence is guaranteed. 

A complete set of invariants for an arbitrary two qubit system as given by Makhlin ll36l is 
listed in Table 12. ll 

It is clear that all the invariants 1^, k = 1, 2, . . . 18, listed in Table 12. ll are invariant under 
local unitary transformations as can be verified by explicitly substituting Eqs. (12.121) . (12.131) 
for transformed state parameters. For example consider the invariant I 4 which, under local 
unitary operation, transform as, 

u = S /T V = s T o« T ow 8 

= s T s since (1)T (9 (1) = 1. (2.14) 

Similarly, it is easy to identify that 



I 12 = s ' T T'r' = s T 0« T O^TO^ T r 
= s T Tr. 



(2.15) 
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h 


= detT 


ho 


= e ijkSi (TT^s) j ([TT' T } 2 S ) k 


h 


= TV (T T T) 


hi 


= e ijhn (T^Tr^dT^T^rh 


h 


= Tr(T T T) 2 


h2 


= s T Tr 


h 


T 

= S S 


hs 


= s T TT T Tr 


h 


= s T TT T s 


hA 


— £ijk £lmn Si T\ tj m t kn 


la 


= s T (TT T ) 2 s 


1 15 


= e- u Si (TT T s)i (Tr)i 


h 


— r T r 


h% 


= (-ijk {T T s)irj (T T Tr) k 


h 


= r T T T T r 


hr 


= tijk (T T s)i (T T TT T s)j r k 


h 


= r T (rr T )2r 


hs 


= e ijk Si (Tr)j (TT T Tr) k 



Table 2.1: Complete set of 18 polynomial invariants for an arbitrary two qubit state. 

Makhlin [36] has given an explicit procedure to find local unitary operations that trans- 
form any equivalent density matrices to a specific canonical form, uniquely determined by 
the set of 18 invariants given in Table |2Tl Further, it has been shown that when T T T is 
nondegenerate, the entire set of invariants Ji_is is required to completely specify the canon- 
ical form of locally equivalent (See Appendix [C]> density matrices. However, when T T T 
is degenerate, only a subset of 18 invariants would suffice for the complete specification 
of density matrices which are locally related to each other. In particular, (i) when two of 
the eigenvalues of T T T are equal then only a subset of nine invariants {I4-9, /12-14} are 
required and (ii) when all the three eigenvalues of T T T are equal, the subset {h-g, I12} 
containing six invariants determines the canonical form of the density matrices. 

In the next section, we identify that the number of invariants required to characterize 
an arbitrary symmetric two-qubit system reduces from 18 (as proposed by Makhlin |[36l ) 
to 6. Moreover, we consider a specific case of symmetric two-qubit system, and show that 
a subset of three independent invariants is sufficient to determine the non-local properties 
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completely. 



2.4 Invariants for symmetric two-qubit states 



Symmetric two qubit states p sym , which obey exchange symmetry, are defined by, 



IIl2 ft™ — ftymni2 — p, 



symj 



where II12 denotes the permutation operator. 

Quantum states of symmetric two qubits get confined to a three dimensional subspace^ 
of the Hilbert space. Explicitly the angular momentum states 

{|f = 1, M = 1), |f = 1, M = 0), |f = 1, M = -1)} are related to the standard two 
qubit states {|0i,0 2 ), |0i, I2), 1 1 1 , 2 ), 1 1 1 , 1 2 )} as follows: 



~ = 1,M = 1 
2 

y = l,M = 



N 



1, M = -1 



|0i, 2 ), 

-^(|0i, l 2 ) + |li,0 2 )), 
|li, h). 



(2.16) 



An arbitrary symmetric two qubit density matrix p sym has the form: 



1 / 3 3 

= - I / % I + Si (an + a 2 i) + Uj o\iOij | , 

i=l ij'=l 



(2.17) 



where, 



Tr(p sym cjH 



Tr(PsymCT2j), 



or r,- 



(2.18) 



'The collective angular momentum basis states {|f , Af); M = — ^ < M < } span the Hilbert space 
"Wsym = (C 2 ® C 2 . . . ® C 2 ) sym of symmetric TV qubit system, i.e., the dimension of the Hilbert space gets 
reduced from 2 N to (TV + 1) for symmetric TV qubit states. 
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and 

Tr[p sym (o-iiO- 2 j)] = Tr[p sym (crijcr 2 i)], 

or Uj = tji T T = T. (2.19) 

(Please see Eqs. (12.61 ) - (12.101) ) for a comparison with the density matrix of an arbitrary two 
qubit system.) 

Obviously, the number of state parameters for a symmetric two qubit system gets re- 
duced from 15 (for an arbitrary two qubit state) to 9 owing to symmetry constraints given 
in Eqs. (12- 18b . (12- 1 9b - Moreover, there is one more constraint on correlation matrix T re- 
ducing the number of parameters required to characterize a symmetric two qubit system to 
8 . To see this, let us consider the collective angular momentum of two qubit system which 
is given by 

J=\(oi+3 2 ). (2.20) 

We have, 

Tr[p sym (J- J)} =2, (2.21) 

for symmetric two qubit states. This is because, the symmetric two qubit system is confined 
to the maximum value of angular momentum j = 1 in the addition of two spin ^ (qubit) 
systems. (Note that in the addition of angular momentum of two spin ^ particles the to- 
tal angular momentum can take values and 1. The two qubit states with net angular 
momentum are antisymmetric under interchange - whereas those with maximum angular 
momentum 1 are symmetric under interchange of particles.) Using Eq. (12.201 ) we obtain, 

J-J = ^[(^1 + <t 2 ) • (<7l + <r 2 )], 

= \[^i-a l ) + {a2-a 2 ) + 2{a l -a2)]. (2.22) 
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Since (<?i • = (0*2 • ct 2 ) = 3/, we can rewrite the above equation as, 



J - J= ^[61 +2 (ffi ■ 3 2 )\. (2.23) 



In other words, we have, 



Tr [p sym (J ■ J) ] = -(3 Tr(p sym J) + Tr[p sym (04 • <t 2 )]) 

= I[3 + Tr(T)]. (2.24) 



From Eq. (12.211 ) we have, 

I[3 + Tr(T)]=2, 
for symmetric two qubit systems. This in turn leads to the constraint 

Tr(T) = l. (2.25) 

Thus only eight real state parameters viz., three real parameters Sj and five parameters Uj, 
(completely specifying the 3x3 real symmetric two qubit correlation matrix T with unit 
trace) determine a symmetric two qubit system. 

We give below an explicit 4x4 matrix form (in the standard two-qubit basis 
I Oi O2) , |0i I2) , I la O2) , 1 1 a I2)) of an arbitrary symmetric two qubit density matrix: 

( l + 2s 3 + i 3 3 A* A* (t u - 1 22 ) - 2it 12 ^ 



1 



A (tll+*22) (tll+*22) B* 

A (t u +t 22 ) (t u +t 22 ) B* 

\ (t n - 1 22 )+ 2it 12 B B l-2s 3 + t 33 J 

(2.26) 

where A = (si + is 2 ) + (ti 3 + it 23 ) and B = (s\ +is 2 )- (h 3 + it 23 ). 
In the following discussion, we show that the number of local invariants required to deter- 
mine a symmetric two qubit density matrix also reduces, owing to the symmetry constraints 
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on the state parameters Eqs. ( 12.181 ), ( 12.191 ). 
Symmetric two qubit local invariants: 

In the case of symmetric two qubit system, it is easy to see that the 1 8 polynomial invariants 
(for an arbitrary two qubit system) given in Table. 12. II in terms of the qubit average s and 
two qubit correlations T reduces to twelve: 



h -- 


= det T 




= h = 


T 

S S, 


h-- 


= Tr(T 2 ) 


h-~ 


= h = 


s T TT T s 


h-- 


= Tr(T 2 ) 2 




= h = 


s T (TT T ) 2 s 




= s T Ts, 




= hi 


= e ijkSi (TT T s),([TT T ] 2 s) k 


h3 


= s T TT T Ts 


hs 


= /l6 


= f-ijkSi (TT T s)j (Ts) k 




^ijk ^Imn Si Si tj m t kn 


hr 


= /l8 


= tijkSi (Ts)j (TT T Ts) k 



Table 2.2: Invariants for a symmetric two qubit state. 

We now proceed to identify that a set containing six local invariants {Zi^q} is sufficient to 
determine the canonical form of locally equivalent symmetric two qubit density matrices. 
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Theorem 2.1 All equally entangled symmetric two-qubit states have identical values for 
the local invariants {Z\ — Iq} given below: 



det 



T, l 2 = Tr{T 2 ) 



I3 = S S , T4 = s 




Z 6 = e ijk Si (Ts)j {T 2 s) k 



(2.27) 



where denotes Levi-Civita symbol; s (s T ) is a column (row) with s\, S2 and S3 as 
elements. 

Proof: Let us first note that the state parameters of a symmetric two-qubit density matrix 
transform under identical local unitary operatiorj^] U (g> U as follows: 

3 

s'j = Oij Sj or s = O s , 



where O € 50(3, R) denotes 3x3 rotation matrix, corresponding uniquely to the 2x2 
unitary matrix U G SU (2). 

To find the minimum number of local invariants required to characterize a symmetric 
two qubit system, we refer to a canonical form of two qubit symmetric density matrix which 
is achieved by identical local unitary transformations U®U such that the correlation matrix 
T is diagonal. This is possible because the real, symmetric correlation matrix T can be 

2 A symmetric state transforms into another symmetric state under identical local unitary transformation on 
both the qubits. 



3 



*ij = 2 OvtOjltu or T' = 0T0 



,T 



(2.28) 



k,l=l 
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It is clear that the invariants 



OTO i 



and To 



h 
t 2 
^ t 3 J 



(2.29) 



2i 
X 2 



det T = t\ t 2 t 3 , 
Tr(T 2 ) = + il + tf, 



(2.30) 



along with the unit trace condition 



Tr(T)=t 1 + t 2 + t 3 = 1, 



(2.31) 



determine the eigenvalues t\, t 2 and t 3 of the two-qubit correlation matrix T. Further, the 
absolute values of the state variables s\, s 2 , s 3 , can be evaluated using T 3 , T4 and X5: 



I3 
X 5 



T 2,2,2 
S s = S 1 + S 2 + S3, 



S T T S = S^ ti + s\ t 2 + S3 t 3 , 

Cijfc Qmn Sj S; tj m tfc n = 2 (sf t 2 t 3 + s\ t\ t 3 + S3 £i £ 2 ) • 



(2.32) 



Having determined s\ , s|, s| and thus fixing the absolute values of the components of the 
qubit orientation vector s - the overall sign of the product S1S2S3 is then fixed by Tq. 



= e ijk Si(TT T s)j ([TT T ] 2 s) k 

= sis 2 s 3 [ti t 2 (t 2 - h) + t 2 t 3 (t 3 - t 2 ) + t 3 ti (ti - t 3 )} . 



(2.33) 



3 Note that Tr(T) is preserved by identical local operations U ®U i.e., we have t\ + ti + ts — 1. 
4 For an arbitrary two qubit state, the diagonal elements (ti, £2, £3) of T d are not the eigenvalues of the 
correlation matrix T (see Appendix [Ct, Hence, to determine (t ly t 2 , £3), (the eigenvalues of TT (TT ) 



local invariants ii_3 are required. However, in the case of a symmetric two qubit state, we have T 
Therefore only two polynomial invariants detT, Tr(T 2 ) suffice to determine the eigenvalues of T. 



T. 
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It is important to realize that only the overall sign of s\ S2 S3 - not the individual signs - is 
a local invariant. More explicitly, if (+, +, +) denote the signs of s±, S2 and S3, identical 
local rotation through an angle it about the axes 1, 2 or 3 affect only the signs, not the mag- 
nitudes of si, S2, S3, leading to the possibilities (+, — , — ), (— , +, — ), (— , — , +). All these 
combinations correspond to the '+' sign for the product si S2 S3. Similarly, the overall '— ' 
sign for the product S1S2S3 arises from the combinations, (— , — , — ), (— , +, +), (+, — , +), 
(+, +, — ), which are all related to each other by 180° local rotations about the 1, 2 or 3 
axes. 

Thus we have shown that every symmetric two qubit density matrix can be transformed 
by identical local unitary transformation U $S> U to a canonical form, specified completely 
by the set of invariants {1i_q}. In other words, symmetric two-qubit states are equally 
entangled iff {Zi^q} are same. In the next section, we consider a special class of symmet- 
ric density matrices and show that a subset of three independent invariants is sufficient to 
characterize the non-local properties completely. 



2.5 Special class of two qubit states 



Many physically interesting cases of symmetric two-qubit states like for e.g., even and odd 
spin states ||49l , Kitagawa - Ueda state generated by one-axis twisting Hamiltonian |[T9l , 
atomic spin squeezed states II201 . exhibit a particularly simple structure 



Qsym 



( a b 

c c 

c c 

y b d 



\ 



(2.34) 



of the density matrix (in the standard two-qubit basis |0i O2) , |0i I2) , |li O2) , |li I2)) with 



Tr(^sym) = a + 2c + d = 1. 



(2.35) 
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The qubit orientation vector s for the density matrices of the form Eq. ( 12.341 ) has the fol- 
lowing structure 

s = (0,0, (a-d)), (2.36) 



V 



and the real symmetric 3x3 correlation matrix T for this special class of density matrix 
has the form: 

2(c + b) 

2{c-b) ■ (2-37) 

(a + d-2c) j 

It will be interesting to analyze the non-local properties of such systems through local in- 
variants. The specific structure £> sym given by equation Eq. ( 12.341 ) of the two-qubit density 
matrix further reduces the number of parameters essential for the problem. Entanglement 
invariants associated with the symmetric two-qubit system g sym , may now be identified 
through a simple calculation to be (see Eqs. (12.30b . (12.321) . (12.331 )) 
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1-2 

2 3 

2 4 

Is 
2 6 



t 1 t 2 t 3 = (Ac 2 -4\b\ 2 )(l-4c), 

4 + 4 + 4 = (2c + 2 \b\) 2 + (2c - 2 |6|) 2 + (1 - 4c) 2 , 

4 + 4 + 4 = (a-d) 2 , 

slh + slt 2 + sjt 3 = (a-d) 2 (l-4c), 

2 (sj t 2 1 3 + s 2 2 h t 3 + S j h t 2 ) = 8 (a - d) 2 (c 2 - |6| 2 ), 

sis 2 s 3 [h t 2 (t 2 - h) + t 2 t 3 (t 3 - t 2 ) + t 3 h (h - t 3 )} = 0. 



(2.38) 



In this special case, we can express the invariants X\ and X 2 in terms of (23,24,15) 
(provided X 3 7^ 0) : 

2 5 2 4 (2 3 -2 4 ) 2 -2 3 2 5 +2| 
11 = Tr2"' 12 = t2 • (2J9 ' ) 

If Z 3 = 0, then the set containing six invariants reduces to the subset of two non-zero 
invariants (2i, T 2 ). Thus the non-local properties of symmetric two-qubit states - having a 
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specific structure g sym given by equation Eq. ( 12.341 ) for the density matrix are characterized 
by 

(i) subset of three invariants (X3, Z4, X5) when X3 7^ or 

(ii) subset of two invariants (Zi,22) when23 = 0. 

In the next section, we propose criteria, which provide a characterization of non- 
separability (entanglement) in symmetric two-qubit states in terms of the local invariants 

{Zi-6>. 

2.6 Characterization of entanglement in symmetric two qubit 
states 

A separable symmetric two-qubit density matrix is an arbitrary convex combination of direct 
product of identical single qubit states, 



(2.40) 



and is given by 

P(sym-scp) = ^2PwPw® Pw, (2-41) 

w 

Where EtB Pv> = l - 

Separable symmetric system is a classically correlated system, which can be prepared 
through classical communications between two parties. A symmetric two qubit state which 
cannot be represented in the form Eq. (12.411 ) is called entangled. 

For a separable symmetric two qubit state, the components of the average spin of the 
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Si — TrfPfsym— sep) 

= Tr[(p w ® p w ) (<rii)] 

w 

= } p w Tr(p w an) 

w 

— ^ ^ Pw Swii (2.42) 

w 

and the elements of the correlation matrix T can be expressed as, 

tij = Tr[p( sym _ sep ) an ® <T2j] 

= ^ Pu, Tr(p w cri i )Tr(p t( ; cr 2 j) 

— ^ Pw S w i S w j- (2.43) 

w 

One of the important goals of quantum information theory has been to identify and char- 
acterize inseparability. We look for such identifying criteria for separability, in terms of 
entanglement invariants, in the following theorem: 

Theorem 2.2 The invariants, Z4, Z5 and a combination Z4 — if of the invariants, neces- 
sarily assume positive values for a symmetric separable two-qubit state, with 1% 7^ 0. 

Proof: (i) The invariant X4 has the following structure for a separable state: 

3 

-^4 — S T S — ^ ^ t"ij Si Sj 

w \i=l J \j=l J 

= (?-^) 2 >Q. (2.44) 

w 
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(ii) Now, consider the invariant X5 for a separable symmetric system: 
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^5 ^ijk tlmn &i S[ tj m tj~ n 

E( (w) (w')\ ( («,) (w') 

Pw Pw' I (-ijk S{ Sj Sfe II f-lmn Si S m S n 



^PwPv 



> 0. (2.45) 



(iii) For the combination X4 — Xf we obtain, 



2 



which has the structure (yl 2 ) — (A) 2 and is therefore, essentially non-negative. Negative 
value assumed by any of the invariants X4, X5 or the X3 — X|, is a signature of pairwise 
entanglement. 

Further from the structure of the invariants in a symmetric separable state, it is clear 
that X3 = ^y^/Pw ys- s^J = implies = for all 'w', leading in turn to X4 = 

w 

ys ■ s^M = and X5 = 0. Thus when X3 = 0, the two qubit local invariant T\ 

w 

characterizes pairwise entanglement as shown in the theorem given below: 

Theorem 2.3 For a symmetric separable two-qubit state, the invariant X\ assumes positive 
value. 

Proof : Consider, 



T = diag(ix, t 2 , t 3 ) 



diag T Pw (s^Y, Tp w (s^)\ T Pw (s^Y) . (2.47) 



We therefore have, 



Xi = det T = h h U = II ( E P» { s i W) )j ' < 2 - 48) 

i=l \ w / 
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which is obviously non-negative for all symmetric separable two qubit states. Thus when 
X3 = 0, X\ > provides a sufficient criteria for separability. 

A simple example illustrating our separability criterion in terms of two qubit local 
invariants is the two qubit bell state 



|$) = _(|01) + |10». 



(2.49) 



The density matrix for the two qubit bell state is written as, 



p = |*)(*| 



0000 
0110 
0110 
0000 



(2.50) 



It is easy to identify that the above density matrix has a structure similar to g sym (see 
Eq. (12.341 )) with the matrix elements given by 



a = 0, 6 = 0, 



c = -, d = 0. 
2 



(2.51) 



The invariant X3 associated with the density matrix of Eq. (12.501 ) is given by, 



2 3 = (a - d) 2 = 0, 



which, in turn implies that 



T 4 = 2 3 (1 - 4c) = 0, 
J 5 = J 3 (c 2 - |6| 2 ) = 0. 



(2.52) 
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T x = (Ac 2 -4|6| 2 ) (1 -4c) = -1. (2.53) 

Since T\ < 0, we can conclude from Theorem 12.31 that the given state Eq. ( 12.491 ) is entan- 
gled. 



2.7 Necessary and sufficient criterion for a class of symmetric 
two qubit states 



It would be interesting to explore how these constraints on the invariants, get related to the 
other well established criteria of entanglement. For two qubits states, it is well known that 
Peres 's PPT (positivity of partial transpose) criterion fl3l is both necessary and sufficient 
for separability. We now proceed to show that in the case of symmetric states, given by 
Eq. (12.341 ), there exists a simple connection between the Peres's PPT criterion and the non- 
separability constraints (see Eqs. (12.441 ) - (12.461 )) on the invariants. 

We may recall from Sec. l2.5l that the density matrix for the special class of symmetric two 
qubit states (see Eq. (12.341 )) has the following structure: 



Qsym 



a b 

c c 

c c 

b d 



\ 



(2.54) 



The partial transpose of the matrix £> sym has the form, 



(o ) PT 



a c 

c b 

b c 

c d 



(2.55) 
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The eigenvalues of the partially transposed density matrix (g S ym) Pr are given by 

Ai = I ((a + d) - V(a - dy + 4 C 2) , 

A 2 = i ((a + d) + V(a - + 4c 2 ) , 

A 3 = c-|6|, 

A 4 = c + |6|, (2.56) 

of which Ai and A3 can assume negative values (a, c, d are positive quantities since they are 
the diagonal elements of the density matrix). 

(i) If Ai < 0, then we have, 

(a + df < {a-df + Ac 2 . (2.57) 
The above inequality can be expressed as 

(a + d + 2c)(a + d- 2c) < {a-df, (2.58) 
which on using Eq. (12.351 ). gets simplified to 

(l-4c) < {a-df. (2.59) 

From Eq. ( 12.381 ), we have 

X 4 - Zf = (a - d) 2 ((1 - 4c) - (a - d) 2 ) , (2.60) 

associated with the two qubit symmetric system g sym . It is obvious that X4 — 2| < 
when (1 — 4c) < (a — d) 2 i.e., when Ai is negative. 

(ii) When A3 < 0, we can easily see that the invariant (see Eq. (12.381 )), 



l 5 = 8(a-d) 2 (c+\b\)(c-\b\) 



(2.61) 
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associated with the special class of symmetric two qubit states (Eq. ( 12.341 )) is also negative. 
Thus A 3 < =>• J 5 < 0. 

We have thus established an equivalence between the Peres 's partial transpose crite- 
rion and the constraints on the invariants for two qubit symmetric state of the form Q( sy m) • 
In other words, the nonseparability conditions I4 — 2|, X5 < 0, are both necessary and 
sufficient for a class of two-qubit symmetric states given by Eq. (12.34b - 
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We have shown that the number of 1 8 invariants as proposed by Makhlin ||36l for an arbi- 
trary two qubit system gets reduced to 12 due to symmetry constraints. Further, we realize 
that a canonical form of two qubit symmetric density matrix achieved by identical local 
unitary transformations U (g> U restricts the minimum number of local invariants to specify 
an arbitrary symmetric two qubit system to 6. In other words, we show that a subset of six 
invariants {Xi_6} of a more general set of 18 invariants proposed by Makhlin ll36l . is suffi- 
cient to characterize the nonlocal properties of a symmetric two qubit states. For a special 
class of two qubit symmetric states, only 3 invariants are sufficient to characterize the sys- 
tem. The invariants X4, T5 and Z4 — Z| of separable symmetric two-qubit states are shown 
to be non-negative. We have proposed sufficient conditions for identifying entanglement 
in symmetric two-qubit states, when the qubits have a non-zero value for the average spin. 
Moreover these conditions on the invariants are shown to be necessary and sufficient for a 
class of symmetric two qubit states. 



Chapter 3 

Collective signatures of entanglement 
in symmetric multiqubit systems 

3.1 Introduction 

Quantum correlated systems of macroscopic atomic ensembles ll50l have been drawing 
considerable attention recently. This is especially in view of their possible applications in 
atomic interferometers Ifl6l [T71 and high precision atomic clocks lfT8l and also in quantum 
information and computation [2]. Spin squeezing [lUCEl has been established as a standard 
collective method to detect entanglement in these multiatom (multiqubit) systems. 

Spin squeezing, is defined as the reduction of quantum fluctuations in one of the spin 
components orthogonal to the mean spin direction below the fundamental noise limit N/4. 
Spin squeezing of around N 10 7 atoms is nowadays routinely achieved in laboratories. It 
has been shown in Ref. [30] that spin squeezing is directly related to pairwise entanglement 
in atomic ensembles - though it provides a sufficient condition for inseparability. 

In the original sense, spin squeezing is defined for multiqubit states belonging to the 
maximum multiplicity subspace of the collective angular momentum operator J. Multiqubit 
states with highest collective angular momentum value J = y exhibit symmetry under 



30 
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the interchange of particles. In other words, the concept of spin squeezing is defined for 
symmetric multiqubit systems and it is reflected through collective variables associated with 
the system. 

In this Chapter, we concentrate on symmetric multiqubit systems and connect the av- 
erage values of the collective first and second order spin observables in terms of the two 
qubit state parameters *] We examine how collective signatures of pairwise entanglement, 
like spin squeezing, manifest themselves via negative values of two qubit local invariants. 
This leads to a classification of pairwise entanglement in symmetric multi-qubit states in 
terms of the associated two qubit local invariants. 

3.2 Collective spin observables in terms of two qubit variables 

The collective spin operator J for a N qubit system is defined by, 

1 - 

J=-J2^a- (3-1) 

Q=l 

Let us concentrate on symmetric multiqubit systems, which respect exchange symmetry: 

n JN) _ JN) jj _ (AO 

where H a p denotes the permutation operator interchanging a th and (3 th qubits. 
The expectation value of collective spin correlations (Jj) is given by 



1 N 

( J i) = 2^^ a ^ ] i = 1 , 2 ,3- 



(3.2) 



a=l 



'As a consequence of exchange symmetry, density matrices characterizing any random pair of qubits drawn 
from a symmetric multiqubit system are all identical. So, the average values of qubit observables are identical 
for any random pair of qubits belonging to a symmetric iV-qubit system. 
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Here, we denote the average value (...) of any observable by 
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(...) = Tr[ Psym (...)]. (3.3) 

It may be noted that the density matrix of the a th qubit extracted from a system of N qubits, 
which respect exchange symmetry, is given by 



(a) rp (AO 

Psym — ir l,2,...a-l,a+l,...7VPsym> 

= \[l + 4 a) Si], (3-4) 



where, 



and 



= J(g>/(g>...(g><7(g>/(g>...(g>/ (3.5) 



is the a th qubit spin operator, with a appearing at a th position. In Eq. (13.41) . the density 
matrix piym of the a th qubit is obtained by tracing the multiqubit state (psym) over all the 
qubit indices, expect a. 

We emphasize that the qubit averages are independent of the qubit index a : 



(^ a) ) = si. (3.6) 



Therefore the first moments of the collective spin (Jj) (see Eq. (13.2b ) assume the form 

1 N N 
(Ji) = 2 = -jSi (3.7) 

a=l 

in terms of the single qubit state parameter Si . 

For our further discussion, we need to associate the second order moments of the 
collective spin variables i.e., (( Jj Jj + JjJi)) with the two qubit con-elation parameters. 
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Using Eq. (13.21 ), we obtain, 



a,P=l 

N 1 N 

= J Sij + - Yl (( a »iCr/3j)) ■ (3.8) 

a^P=l 

Here, {{cr a i&pj)) (a 7^ /?) are the spin correlations of a pair of qubits a, j3 drawn randomly 
from a symmetric multiqubit system. 

The density matrix (psym ) of sucn a P arr of qubits obtained by taking a partial trace 
over the remaining (N — 2) qubits is given by 

(Psym ) Tri,2, .. .except (a,/3) (Psym) - 

The general form of such a two qubit density matrix in terms of 8 (see Eq. (12.171) ) state 
variables is given by, 



/ N N 

pit = \ I 1 ® 1 + E « + ^) + E ^ °i 0) ^ I 

\ i=l a,P=l 



where, 



= <(^Vf)>, (3.10) 



irrespective of the qubit indices a, /?. 

Substituting Eq. (13.101 ) in Eq. (13.81 ). we obtain 



1 N N(N-l) 

-((JiJj + JjJi)) = + -((W2j)) 



N 

■j[Si3 + (N-l)Uj], i,3 = 1,2,3. (3.11) 
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From Eqs. ( 13.71 ), (13.111 ), it is evident that the collective spin observables (upto first and sec- 
ond order) can be expressed in terms of state parameters of a pair of qubits chosen arbitrarily 
from a symmetric N qubit state. Thus the collective pairwise entanglement behavior in sym- 
metric multiqubits results from the properties of the two qubit state parameters {sj, tij} . 

We now proceed to identify collective criteria of pairwise entanglement in a symmetric 
multi-qubit state, in terms of the two-qubit local invariants {T\ — Xq}. 

3.3 Collective signatures of pairwise entanglement 

Collective phenomena, reflecting pairwise entanglement of qubits, can be expressed through 
two qubit local invariants as the first and second moments ( Jj), ((JiJj + JjJi)) are related 
to two qubit state parameters {sj, t^} in symmetric multiqubit systems. Here, we show 
that spin squeezing- which is one of the collective signatures of pairwise entanglement 
in symmetric multiqubit systems-gets reflected through one of the separability criterion 
derived in Chapter |2 

For the sake of continuity, the main result of Chapter|2]is summarized in the following 
sentence: 

The non-positive values of the invariants {1$, I§ or Z4 — j|} serve as a signature of 
entanglement and hence provide sufficient conditions for non-separability of the quantum 
state. 

Let us now review spin squeezing criteria. Kitagawa and Ueda [19 ] pointed out that a 
definition of spin squeezing, |5T| based only on the uncertainty relation, 



(AJ!) 2 (AJ 2 ) 2 > 



(3.12) 



4 



exhibits co-ordinate frame dependence and does not arise from the quantum correlations 
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among the elementary spins. They identified a mean spin direction 

* = (3.13) 



where, \ (J)\ = y (J) • (J) (The collective spin operator J for an N qubit system is given 
by Eq. (IPT ). 

Associating a mutually orthonormal set {h\±, n 2 j_, n } , with the system, let us con- 
sider the following collective operators, 

Ji_l = J ■ ni±, J 2 j_ = J • n 2 j_ and Jq = J ■ uq (3.14) 

which satisfy the usual angular momentum commutation relations 

[Ju_,J 2 ±]=iJo. (3-15) 

Now, employing a collective spin component J± orthogonal to the mean spin direction no, 
given by, 

J_l_ = J • nj_ 

= J\± cos 9 + J 2 j_ sin 6, (3. 16) 

minimization of the variance, 

(AJ ± ) 2 = (Jl) - (J ± ) 2 (3.17) 



can be done over the angle 9. Kitagawa and Ueda Ifl9ll proposed that a multiqubit state can 
be regarded as spin squeezed if the minimum of A J± is smaller than the standard quantum 
limit of the spin coherent state. 
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A spin squeezing parameter incorporating this feature is denned by |[T9l 



2(AJ ± ) 



(3.18) 



Symmetric multiqubit states with £ < 1 are spin squeezed. We next proceed to show that 
the two qubit local invariant J5 and the spin squeezing parameter £ are related to each other. 

3.3.1 Spin squeezing in terms of the local invariant J 5 

We now prove the following theorem. 

Theorem 3.1 For all spin squeezed states, the local invariant Z5 is negative. 

Proof. It is useful to evaluate the invariant X 5 (see Eq. ( 12.271 )), after subjecting the quantum 
state to a identical local rotation U ®U (8> U (g> • • • on all the qubits, which is designed to 
align the average spin vector (J) along the 3-axis. After this local rotation, orientation of 
the qubits would be along 3 axis and the qubit orientation vector is given by 



s = (0, 0, s ). 



We may then express the local invariant X5 as, 



X5 




■n 



^3jk Qmn tj m t)~ n 
2 Sq (t 11*22 — (t'u) 2 ) 



2 det T 



(3.19) 



where, 



(3.20) 
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denotes the 2x2 block of the correlation matrix in the subspace orthogonal to the qubit 
orientation direction i.e., 3-axis. 

Now, we can still exploit the freedom of local rotations O12 in the 1 — 2 plane, which 
leaves the average spin s = (0, 0, sq) unaffected. We use this to diagonalize T±: 



12 T l OTo = T 



( 



12 



(3.21) 



with the diagonal elements given by 



t (±) 



(t' n + t' 22 )±J(t' n -t' 22 y + A(t> 12 y 



We once again emphasize that local rotations on the qubits leave the invariants unaltered 
and here, we choose local operations to transform the two qubit state variables as, 



(0,0, s ), 



(3.22) 



and 



( 4 +) 



13 



1 23 



(3.23) 



\ l 13 z 23 l 33 ) 



so that the two qubit invariant X5 can be expressed as, 



T 5 = 2s§detTj_, 
_ 2 2 ,(+),(-) 



(3.24) 



Collective signatures of entanglement in symmetric multiqubit systems 



38 



We now express the spin squeezing parameter £, given by Eq. ( 13.181 ), in terms of the 
two-qubit state parameters 



4(AJ ± ) 



2 

min 



N 

4({J-fi±) 2 )-((J-fi±)Y' 
N 



(3.25) 



Writing the collective spin operator J for an N qubit system in terms of the two qubit 
operators (see Eq. (13.11 )). we obtain 



1 N 

= N Yl ((°a-h ± )(a f 3-h ± )) i 



a,/3=l 



N N 



1 + J^J2 12 ((^■n±)(a p -h±)) I 



a=l f3^a=l 
N N 3 



a=l f3>a=l \i,j=l 



(3.26) 



Since for a symmetric system, we have (cr a i&pj) = Uj, we express the squeezing parameter 
as follows: 



£ 2 = 1 + (JV - 1) \J2 Uj n ±i n ±j 
= 1 + (N- l)(nlTn ± ) min . 



(3.27) 



In Eq. (I3.27I ). we have denoted the row vector = (rii±, ri2±, 0) = (cos 9, sin#, 0) . 
The minimum value of the quadratic form (n^ Tn_|_) m i n in Eq. (13.271 ) is fixed as follows: 



( ■ \ 

mm 



(t'n cos 2 9 + t' 2 2 sin 2 9 + t' 12 sin 29) 



— (t'n + ^22) 



(t il -t 22 ) 2 + 4(V 2 



(3.28) 
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where t\_ is the least eigenvalue of T± (see Eq. (13.201 )). 



We finally obtain, 



£ 2 = ±(AJ ± ) mi „= (i + (iv-i)t; 



(-) 



) 



(3.29) 



Following similar lines we can also show that 




max 



l + (iV-l)^ +) ), 



(3.30) 



which relates the eigenvalue t ± of T± to the maximum collective fluctuation (A Jj_) ma: 
orthogonal to the mean spin direction. Substituting Eqs. d3.291 ), (13.30I ). and expressing 



Having related the local invariant T5 to collective spin observables, we now proceed to show 
that X5 < iff £ 2 < 1 i.e., iff the state is spin squeezed. 

Note that the two qubit correlation parameters Uj are bound by 



This bound, together with the unit trace condition Tr (T) = 1 on the correlation matrix of a 
symmetric two-qubit state, leads to the identification that only one of the diagonal elements 
of T can be negative. This in turn implies that if the diagonal element t\_ ' is negative, then 
the other diagonal element is necessarily positive. Thus, from Eq. ( 13.301 ), it is evident 
that 



so = jf\(Js) I in Eq. {HQ), we get, 




(3.31) 



-1 < Uj < 1- 
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whenever < 0. It is therefore clear (from Eq. (13.311 )) that a symmetric multiqubit state 
is spin-squeezed iffls < 0. In other words, 

£ 2 < 1 <^=> X 5 < 0. (3.32) 

Further, from the structure of the invariant X5 (Eq. ( 12.271 )), it is clear that X5 < implies 

i.e., one of the eigenvalues t\, ti or ts of the correlation matrix T must be negative. This in 
turn implies that the invariant 

Ji = ht 2 t 3 < 0. 

In other words, when X5 < 0, the invariant T\ is also negative. 

We now explore other collective signatures of pairwise entanglement, which are man- 
ifestations of negative values of the invariants X4 and X4 — X| . 



3.3.2 Collective signature in terms of J 4 

When the average spin is aligned along the 3-axis through local rotations such that 
s = (0, 0, so)> the local invariant X4 (see Eq. (12.271 )) assumes the form, 



Ia = s T s 



s 



/ +(+) 
z ± 



t" 







t" \ 

c 13 



t { ~ ] t" 



t" 

E 23 



33 / 





V 50 J 



s 2 Q t' 33 . (3.33) 



It is evident from the above equation that 



X 4 < iff tL < 0. 
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Simplifying Eq. ( 13.71 ) and Eq. ( 13.111 ), we express t 33 and so in terms of the collective spin 
observables i.e., 

2 _ 2 

4 = 7^1) (|((^-o) 2 )-l), (3-34) 
leading further to the following structure for the invariant Zj 

T * = NHN-1) l(/)|2 {j?« f - ^ ~ l ) ' (3 " 35) 

where no denotes a unit vector along the direction of mean spin. We therefore read from 
Eq. (13.351 ). that the average of the squared spin component, along the mean spin direction , 
reduced below the value Nf 4, signifies pairwise entanglement in symmetric N-qubit system. 

Further, we may note that when Z4 < 0, the invariant Z5, which reflects spin squeezing 
is not negative . This is because t 33 < => t± > as, 

(i) + tj_ + t' 33 = 1, (unit trace condition), and 

(ii) -1 t 33 < !■ 

Therefore, spin squeezing and ((J ■ no) 2 ) < -j are two mutually exclusive criteria of 
pairwise entanglement. 

However, from the structure of the invariant Z4, as given in Eq. (12.321) . it is obvious 

that 

1 A = slh + s 2 2 t 2 + s 2 3 t 3 <0 (3.36) 

implying that one of the eigenvalues t\, t 2 , £3 of two qubit correlation matrix must be neg- 
ative. This leads to the identification that 

Zi = h h h < 0. 



We now continue to relate the combination of invariants X 4 — X| to the collective spin 
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observables. 



3.3.3 J 4 — If in terms of the collective variables 



In terms of the two qubit state parameters, expressed in a conveniently chosen local coor- 
dinate system (see Eqs. ( 13.221 ), ( 13.231 )), the invariant combination X 4 — Xf can be written 
as, 

X 4 - Xf = s T T s - s T s = s 2 (4 - ag). (3.37) 

Since the state parameters can be expressed in terms of the expectation values of the collec- 
tive spin variables Eq. (13.34b . the invariant quantity X 4 — X| may be rewritten as 



I A — Xo 



N 2 



N(N - 1) 



((J -n ) 2 



1 



(N - 1) N 2 



16 



N 3 (N - 1) 



\(J)\'' 



((J -n ) 



N ' (N - 1] \(J)\ 2 



4 + 



N 



(3.38) 



Negative value of the combination X4 — X| manifests itself through 



From Eqs. (13.351 ) and (13.381 ). we conclude that pairwise entanglement resulting from 



X 3 / 0, X 4 > 0, but X 4 - Xf < 0, 



realized, whenever 

- < ((J ■ n ) 2 ) < - + - |(J)| 2 . 

All the cases discussed above are valid when the average spin vector \ (J)\ / i.e., when 
I (j) I is oriented along the 3-axis. In the special case when | ( J) \ = 0, we show that pairwise 
entanglement manifests itself through negative value of the local invariant X 4 . 
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3.3.4 Characterization of pairwise entanglement through Z x 

In the cases where the qubits have no preferred orientation, i.e., when | (J) | = 0, it is evident 
from Eq. ( 12.271 ) that the local invariants {Tz_q} are zero 

J 3 = S T s = 0, 
1 4 = s T Ts = 0, 

2-5 — tijk tlmn Si Si tj m t kn — 0, 

J 6 = e ijk s i (Ts) j (T 2 s) k = 0. (3.39) 
The remaining two nonzero invariants T\ and Z 2 are, 

X\ = det T, 

2" 2 = Tr(T 2 ). (3.40) 

In such situations, i.e., when \(J)\ = 0, pairwise entanglement manifests itself through 
Ji < 0. 

Writing the invariant X\ in terms of collective observables Eq. (13.111 ). we have, 

1, = detT = tl t 2 t 3 = (^TT)) 3 ft (<J?> " t) • ( 3 - 41 ) 

Negative value ofT\ shows up through (J?) < ^ along the axes i = 1, 2 or 3, which are 
fixed by verifying ((JjJj + JjJi)) = 0; i / j, as T is diagonal with such a choice of the 
axes. 

Note that, 

Ji = det T < J 2 = Tr (T 2 ) > 1 

since t\ + t 2 + *3 = 1 (unit trace condition) and — 1 < t\, i 2 , *3 < 1- Therefore we have, 



X! < 0, J 2 > 1, 



(3.42) 
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both implying pairwise entanglement. 

Thus, we have related the two qubit entanglement invariants Eq. (12.27b to the collective spin 
observables and shown that the collective signatures of pairwise entanglement are mani- 
fested through the negative values of the invariants Z4, Z5, Z4 — X|. 
In Sec. 13.41 (a) we relate the invariant criteria with the recently proposed generalized spin 
squeezing inequalities for two qubits 0T1 and (b) propose a classification scheme for pair- 
wise entanglement in symmetric multiqubit systems. 



3.4 Classification of pairwise entanglement 

Recently, Korbicz et al. 11311 proposed generalized spin squeezing inequalities for pairwise 
entanglement, which provide necessary and sufficient conditions for genuine 2-, or 3- qubit 
entanglement for symmetric states: These generalized spin squeezing inequalities are given 
by EH 

W<l-« (343) 
where Jy. = J • k; with k denoting an arbitrary unit vector. 

We now show that that the generalized spin squeezing inequality given in Eq. ( 13.431 ) can be 
related to our invariant criteria. 

We consider various situations as discussed below: 
(i) Let k = nj_, a direction orthogonal to the mean spin vector (J). The inequality given by 
Eq. (13.431 ) reduces to 

(AJ„J 2 < j. 
Minimizing the variance (AJ n± ) 2 we obtain 

2 4(AJ ± )^ n 
^ = N <L 
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This is nothing but the conventional spin squeezing condition T5 < in terms of the invari- 
ant. 

(ii) If k is aligned along the mean spin direction i.e., k = h$ with no = 77% > the generalized 
spin squeezing inequalities (see Eq. ( 13.431 )) reduce to the form, 

((/•n ) 2 )<f + ^^|(J)| 2 . (3-44) 



From Eq. (13.351 ) we have, 



\{J)\ 2 ( 4<(^o) 2 > - 1 ) • (3.45) 



* N 2 {N- 1) IN 71 

Now the condition T4 < on the local invariant leads to the collective signature [see Ta- 
bid] 

<(J ■ no) 2 } < T , 



which is a stronger restriction than that given by Eq. (13.44b . 

3 



Further, if T 4 > but X 4 - Jf < we obtain the inequality [see Table. 1] 



- < ((J ■ n ) 2 ) < j + - J \(J)\ 2 , 



which covers the remaining range of possibilities contained in the generalized spin squeez- 
ing inequalities of Eq. (13.441 ) with n along the mean spin direction. 

(iii) If the average spin is zero for a given state i.e., we have ( = for all directions k. 
The inequalities of Korbicz et al. OTTl assume a simple form 



« 2 < t 



This case obviously corresponds to X3 = and X\ < 0. 
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In the following table we summarize the results and prescribe a classification of pair- 
wise entanglement in symmetric multiqubit states. 



Criterion of pairwise entanglement 


Collective behaviour to look for 




J 5 <o 


(AJj_)^ in < f 




J 4 < 






1 4 > 0, 1 4 - If < 


f«(Mo) 2 )<f + ¥l(f 


T 3 = 


li < 


(J?) < f 
for any direction i = 1, 2, 3, so that 

((Ji Jj + Jj Ji)) = 0; for i + j 



Table 3.1: Classification of pairwise entanglement in symmetric multi-qubit states in terms 
of two-qubit local invariants. 
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In summary, we have shown that a set of six local invariants {T\ — Iq}, associated with the 
two-qubit partition of a symmetric multiqubit system, characterizes the pairwise entangle- 
ment properties of the collective state. We have proposed a detailed classification scheme, 
for pairwise entanglement in symmetric multiqubit system, based on negative values of the 
invariants X±, Z4, X5 and X4 — X| . Specifically, we have shown, collective spin squeezing 
in symmetric multi-qubit states is a manifestation of X5 < 0. Moreover, we have related 
our criteria, which are essentially given in terms of invariants of the quantum state, to the 
recently proposed generalized spin squeezing inequalities OH for two qubit entanglement. 



Chapter 4 



Dynamical models 



In the previous chapters, we have proposed separability criteria for symmetric multiqubit 
states in terms of two qubit local invariants. In the light of our characterization for pairwise 
entanglement, we analyze few symmetric multi-qubit dynamical models like, 

1. Dicke states 0911401 

2. Kitagawa-Ueda state generated by one axis twisting Hamiltonian |fT9l 

3. Atomic squeezed states |[20l . 

4.1 Dicke State 

Collective spontaneous emission from dense atomic systems has been of interest since the 
pioneering work of Dicke, who predicted that two-level atoms (or qubits) possess collective 
quantum states in which spontaneous emission is enhanced (superradiance) or suppressed 
(subradiance). Multiqubit Dicke states are of interest for quantum information processing 
because they are robust under qubit loss ll52l l53l and stand as an example of decoherence- 
free subspaces. An N-qubit symmetric Dicke states | J = y,M) ; —J<M<J, with 
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(y — M) excitations (spin up) is denned as ll52l 



N 



) 



/ u Pk (| l-i 3 12, ••■1m, m+ i, ...Oa?)) 



(4.1) 



V 



M 



A' 



where m = y — M and {-Pfc} is the set of all distinct permutations of the spins. A well 
known example is the W-state given by, 



— , M = — - 1 ) = -= [|li0 2 3 ■ • ■ Ojv> + IO1I2O3 • • • Ojv) + ■ • • + |0i0 2 3 • • • l N )] ,(4-2) 



which is a collective spin state with one excitation. 

Multiqubit Dicke states are symmetric under permutation of atoms and entanglement- 
robust against particle loss f52~1 . They exhibit unique entanglement properties ll54l and are 
excellent candidates for experimental manipulation and characterization of genuine multi- 
partite entanglement. It has been shown that a wide family of Dicke states can be generated 
in an ion chain by single global laser pulses ll55l . Further, a selective technique that allow 
a collective manipulation of the ionic degrees of freedom inside the symmetric Dicke sub- 
space has been proposed ll56l . An experimental scheme to reconstruct the spin-excitation 
number distribution of the collective spin states i.e., tomographic reconstruction of the diag- 
onal elements of the density matrix in the Dicke basis of macroscopic ensembles containing 
atoms, with low mean spin excitations, has also been put forth ll57l . More recently, Thiel et. 
al. ll58l proposed conditional detection of photons in a Lambda system, as a way to produce 
symmetric Dicke states. 

In order to analyze the pairwise entanglement properties of N-qubit Dicke states in 
terms of two qubit local invariants, we need to evaluate the first and second order moments 
(Jj), (JiJj + JjJi) of the collective spin observable. These moments in turn allow us to 
determine the two qubit state parameters associated with a random pair of qubits (atoms), 
drawn from a multiqubit Dicke state. 



N 



2 



N 



2 
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Average values of the collective spin observable J: 

It is easy to see that 

(J 1 ) = (J,M\J l \J,M) = ±[(J,M\J + + J_\J,M)} 

= 

(J 2 ) = (J,M\J 2 \J,M) = ±[(J,M\J+-J-\J,M)] 

li 

= 

(J 3 ) = (J,M\J 3 \J,M) = M[(J,M\J,M)] 

= M. (4.3) 

The average values of second order collective spin correlations in N-qubit Dicke states are 
readily evaluated and are given by, 

(J 2 ) = {J,M\jf\J,M) 

= ^(J,M\(J + + J^) 2 \J,M) 

= - \(N 2 + 2N - AM 2 )] 
8 

(N 2 + 2N - 4M 2 ) 



(Jf) = (J,M\J 2 \J,M) 

= ±(J,M\(J + -J^) 2 \J,M) 

= - \(N 2 + 2N - 4M 2 )] 
8 

(iV 2 + 2N - AM 2 ) 
8 ' 



= (J,M\J 2 \J,M) 
= M 2 (J,M\J,M) 
= M 2 , 



(4.4) 
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(Ji J 2 + J 2 Ji) 



1 

T 
o. 



■(J,M\(J + + J_)(J + - J_) + (J+ + J_)(J 4 



51 

J_)|J,M) 

(4.5) 



Similarly, we find that 



(JiJj + JjJi) = for i ± j. 



(4.6) 



4.1.1 Two qubit state parameters for Dicke state 

We may recall here that the components of the single qubit orientation vector Si (Eq. (12.91 )) 
are related to the first order moments of the collective spin observables (see Eq. ( 13.71 )) 
through the relation (Ji) = ySj. Thus it is clear from Eq. (14.31) . that 



51 = iv (Jl) = 0, 



S2 = jy(J2>=0, 



2 2M 

53 = iv (Js) = l^- 



(4.7) 



In other words, the qubit orientation vector of any random qubit drawn from a iV qubit 
Dicke state has the form, 



2M 

^(0, 0,— 



(4.8) 



As the qubit correlations are related to the collective second moments (JjJj + JjJi 
through (see Eq. (13.111 )) 



tij ~ N - 1 



2{J i J j + J j J i ) 
N 



'i j 



(4.9) 
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we can evaluate the matrix elements of T by using Eq. ( 14.4 



4[(J! 2 )] 1 
N(N-l) N-l 



1 



4{N 2 + 2N - AM 2 



N(N - 1) 
N 2 - AM 2 
2N(N - 1)' 



N 



t22 



4[(Jf)] 1 
N(N-1) N-l 



1 



4(7V 2 + 2iV - 4M 2 



iV(iV - 1 
N 2 - AM 2 
2N(N - 1) ' 



2V 



*33 



1 



4[(Jp] 
iV(iV - 1) AT - 1 

4M 2 - 2V 
iV(iV - 1)' 



(4.10) 



Further from Eq. ( 14.61 ), it can be easily seen that the off diagonal elements of the correlation 
matrix corresponding to Dicke state are all zero i.e., 



tij = with i ^ j. 



(4.11) 



We thus find that the 3x3 real symmetric two qubit correlation matrix T associated for any 
random pair of qubits, drawn from a multiqubit Dicke state is explicitly given by, 



/ 



T = diag (ti, t 2 , t 3 ) 



N 2 -4M 2 
2N(N-1) 








N 2 -4M 2 
2N(N-1) 











\ 



4M 2 -N 
N(N-l) / 



(4.12) 



Dynamical models 



53 



We construct the density matrix characterizing a pair of qubits arbitrarily chosen from a 
multiqubit Dicke state 



Qsym 



' a 

c c 

c c 

y o o o d 



\ 



where, 



(N + 2M)(N -2 + 2M) 



N 2 



4N(N 
4M 2 



1) 



4N(N - 1) ' 
(N - 2M)(N 



2M) 



4N(N - 1) 



(4.13) 



(4.14) 



We may note here that the above density matrix belongs to the special class of density 
matrices discussed in Sec |2.5i Thus the non-local properties of symmetric N-qubit Dicke 
states are characterized either by subset of three invariants (T^,T^,T^), when T3 / or a 
subset of two invariants (Tx,^), when 2" 3 = 0. We now proceed to evaluate the two qubit 
local invariants associated with the two qubit density matrix Eq. ( 14. 131 ). 



4.1.2 Local invariants 

The two-qubit local invariants (Eq. ( 12.271 )), associated with Dicke state can be readily eval- 
uated and we obtain, 



Ti = det T = ti t 2 h 

f N 2 - AM 2 \ 2 f AM 2 -N\ 
\2N(N -1) ) \N(N - 1)/ ' 



(4.15) 
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Tr (T 2 ) = t\ + t\ + t§ 



/ iV 2 - 4M 2 \ 2 / 4M 2 - N 



2 \2N(N -I)) + \N(N-l] 



T 2,2,2 
S S = S 1 + S 2 + S3 



4M 2 



S T T S = S 2 t\ + S 2 *2 + «3 *3 



( AM 2 - N 
\N(N - 1) 



Cij'fc f-lmn Si Si tj m ifc n 
2(sl t 2 t 3 + s|tlt 3 + s|tlt2) 
8J3 



TV 2 - 4M 2 x 2 



4N(N-1)J ' 



e iife Sj (Ts)j (T 2 s) fc 

SiS 2 S 3 [*i *2 (*2 - *i) + h t 3 (t 3 - t 2 ) + t 3 ti (tx - t 3 )] 

0. (4.16) 



Further, the combination T4 — I 2 of invariants, is given by 

9 ( AM 2 — N 2 \ 

We now consider three different cases (for different values of M) and explicitly verify 
pairwise entanglement of Dicke states through two qubit local invariants, 
(i) When M = ±f : 

In this case, the multiqubit dicke state state has the form 

N N\ 

y, y) = |0i0 2 • • -Oat) . (4.18) 
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This corresponds to a situation in which all the qubits are spin-up . The N-qubit Dicke state 
in which all the qubits are spin-down is given by. 



N N\ 

y ~~2/ 



iii 



]±2 • • • 



In) 



(4.19) 



The collective state, corresponding to this case, is obviously a uncorrected product state. 
The invariants in this case are given by 



T\ — T5 
I2 =X 3 



0. 

^ = 1, 



(4.20) 



which are all non-negative indicating that 1^, Dicke states are separable, 
(ii) M=0: 



TV 
2 ' 



M = 0) Dicke states are written as 



N 



lil 



1-1-2 • 



1jv,0jv_|_^0jv_|_2 ' " "O^v 



(4.21) 



The invariants in this case are given by 



J 3 = J 4 = X 5 = 0, 



(4.22) 



while the non-zero invariant, 



I, 



1 



N 



4 V N- 1 



(4.23) 



assumes negative value. 



So, the Dicke state 1^, 0), (with even number of atoms), exhibits pairwise entanglement, 



which is signalled in terms of collective signature (see Table. (13.11) ) (Jf) < 



N 
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(iii)M^if , 0: 



2 

In this case, the invariant, Z4 is bound by 



N-l 



<1 4 <1, 



and the combination Z4 — 2| is always negative, thus revealing pairwise entanglement in 
Dicke atoms in this case too. The corresponding collective signature is given by -j < 
((J - no?) < f + ^ I (J) | 2 (see Table. (ED). 



4.2 Kitagawa-Ueda state generated by one axis twisting Hamil- 
tonian 

In 1993, Kitagawa and Ueda flT9l had proposed the generation of correlated iV-qubit states, 
which are spin squeezed, through the nonlinear Hamiltonian evolution H = j\ x> 

|*K-u) = e~ lHt I J, - J) ; J = y , (4.24) 

referred to as one-axis twisting mechanism. The N qubit state | J, — J) is the all spin down 
state 

\J,—J) = 1 1 1 , 12, 13, •■•Ijv)- 

The one-axis twisting Hamiltonian has been realized in various quantum systems in- 
cluding quantum optical systems |[59l . ion traps |[60l . cavity quantum electro magnetic dy- 
namics lloTll . This effective Hamiltonian H = j\ x> nas already been employed to produce 
entangled states of four qubit maximally entangled states in an ion trap [62). Collisional in- 
teractions between atoms in two-component Bose-Einstein condensation are also modeled 
using this one-axis twisting Hamiltonian lTT3l . 

In order to investigate the entanglement properties for a random pair of qubits drawn 
from the Kitagawa-Ueda state, we evaluate the first and second order moments of the 
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collective spin operators. 

Expectation values of the spin operator J: 

To evaluate the expectation values (Jj) and (JiJj + JjJi), let us first consider the time 
dependent operators J 3 (t) under the Hamiltonian evolution: 

J 3 (t) = e iHt he~ lHt = J 3 + [iHt, J 3 ] + ^[iHt, [iHt, J 3 ]] + ... (4.25) 

(Here we have used the Baker-Campbell-Hausdorf formula e A Be~ A = B + [A, B] + 
1,[A,[A,B]} + ^[A, [A, [A, £]]] + •••) 

The commutators in Eq. (14.251) are given by: 

[iHt, J 3 ] = [i X tJl J 3 ] = ixt {Ji [Ji, Js] + [Ji, Js] Ji} 
= ixt{Ji(-iJ2) + {-iJijJi} 

where we have denoted [0\, C*2]+ = O1O2 + OiO\. We further obtain, 

[iHt, [iHt, J 3 ]] = [ix^i,X^i,</ 2 ]+] 

= i X 2 i 2 {Ji[Ji, [Ji, J 2 ]+] + [Ji, [Ji, J 2 ]+] <M 
= -x 2 * 2 { Ji J3 + 2 Ji J 3 Ji + J 3 J 2 } • 

Therefore we get, 

■h (t) = J 3 + Xt[Ji , J 2 ]+ - ^X 2 * 2 { Ji 2 ^3 + 2 Ji J 3 Ji + J 3 J x 2 } + ■ ■ ■ • (4-26) 
Now we consider the time dependent operator J+(t) under the Hamiltonian evolution H, 

J+ (t) = e iHt J + e~ iHt = J+ + [iHt, J+] + - [iHt, [iHt, ■/+]] + •••. (4.27) 
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The commutators can be computed as follows: 

[iHt, J+] = [i X tJl J+] = i X t {Ji[Ji, J+] + [Ji, J+] Ji} 
= »x*{Ji(-J 3 ) + {-Js)Ji} 
= -ixt[J\,h] + 

[iHt, [iHt, J + }] = [ixtJl-ixt[Ji,J 3 }+} 

= X 2 t 2 {Ji[Ji, [Ji, Js]+] + [Ji, [Ji, J 3 ]+] Ji} 
= -i X 2 t 2 {Jlh + 2 Ji J 2 Ji + J 2 Jf} . 

Hence, 

J+(t) = J+- i X t[Ji, J 3 ]+ - iX 2 t 2 {Jlh + 2 hhh + J 2 Jf} + -- - ■ (4.28) 
Similarly J_ (t) can be evaluated as 

J_(t) = e iHt J_ e - iHt = J_ + [iHt, J_] + ^[i-ffi, [i£Ti, J_]] + • • • (4.29) 

[itft, J_ ] = [i X tJlJ-]=ixt{ Ji [Ji , J-] + [Ji, J-]Ji} 
= ixt{Ji{-h) + {-h)Ji) 
= ixt[Ji,M + 



[iHt, [iHt, J-]] 



= [ixtJhixt[ J i^MA 

-x 2 t 2 {Ji[Ji, [Ji,J 3 ]+] + [Ji,[Ji, J 3 ]+]Ji} 
= i X 2 t 2 { Jl J 2 + 2 Ji J 2 Ji + j 2 jI) . 
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So we obtain, 



J_ (t) = J_ + ix*[Ji, J 3 ]+ + «X 2 t 2 {J? ^2 + 2 Ji J 2 Ji + J 2 J? } + ■■■ ■ (4.30) 



Using Eqs. (14.26b . (14.281 ). (14.30b . the expectation values are evaluated as follows: 



<*K-u|./3(t)|*K-U> = (J ~ J\ J 3 (t)\J ~ J) 



1 



(J ~ A [Js + Xt[Ji, J 2 ]+ - ^X 2 t 2 { Jlh + 2 Ji J 3 Ji + JzJi) + • • • ] I J - J) 

2V 



21 

1 ,.„„ _ 1 



2 

iV 
— cos 
2 



l-^X 2 t 2 {N-l)] + - X H\N-lf 



N-l 



(xt), 



(4.31) 



(*K-u|^+(t)|*K- 



{J - J\J+(t)\J - J) 

(J -J\J+- i X t[Ji, J 3 ]+ - ix 2 t 2 {J?J2 + 2 Ji J 2 Ji + J 2 J x 2 } + • • • | J - J) 
(4.32) 



<* K -u|J-(t)|*K-U> 



(J — J|J_(t)| J — J) 

(J - J| J_ + *x*[Ji, J 3 ]+ + ix 2 * 2 { J?^ + 2 Ji J 2 Ji + J-2 J 2 } + • • • | J - J) 
(4.33) 



(^k-u|«/ 3 2 |^k-u) 



(J — J| jf (t)|J — J) 
(J- J| ( J 3 + xt[Ji,J 2 ] 



-[iV 2 + iV + iV(iV - 1) cos JV - i (2xt)] 



^yX 2 i 2 {JfJa + 2 Ji J 3 h + J 3 J 2 } + 



| J- J) 
(4.34) 



In a Similar manner, the remaining first and second order expectation values can be ob- 
tained. The average values of the collective spin observables for Kitagawa-Ueda state are 
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listed below: 



<Jl> = (J 2 ) = o, 

(Js) = -ycos^xt) 



(J!) 



N 

T' 

(J|) = i(iV 2 + iV-iV(iV-l)cos Ar " 2 (2 X t)) 
(J|) = I(Ar2 + iV + i v (i v_i) cos ^2 (2x ^ ; 

<[Jl,J 2 ]+) = ^iV(iV- l)^- 2 (xt)sin(xt), 

<[J + ,J 3 ] + ) = 0. (4.35) 



4.2.1 Two qubit state variables 

The qubit state parameter Si associated with a random qubit chosen from a multiqubit 
Kitagawa-Ueda state can be written (from Eq. ( 14.351 )) as, 

si = |fW = o, 

S3 = ^(J 3 ) = -co^ N - l \ X t). (4.36) 

Therefore the orientation vector s for qubit drawn randomly from the Kitagawa-Ueda state 
is given by 

s=(0, 0, -cos^-^xt)). 
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The two qubit correlation matrix elements are readily obtained from the second order mo- 
ments of the collective spin observable Eq. (14.351) . 



1 



N - 1 



2 i^JiJj -\- JjJi 
N 



Thus we have, 



hi 

t-22 
h3 



tl3 — *23 — 0, 
cos^-^xi) sin(xt), 
\ (l-cos^ 2 )(2 X t) 
\ (l + cos(^ 2 )(2 X t) 



(4.37) 



The 3x3 real two qubit correlation matrix T can be thus written as, 



\ 



V 



co$( N -V(xt) sin(xt) 

cos^- 2 ) (x t) sm( X t) \ (1 - cos^" 2 ) {2 X t)) 

i (l + cos(^ 2 )(2 X t)) ) 



The two qubit density matrix for a pair of qubits arbitrarily chosen from a multiqubit 
Kitagawa-Ueda state is given by, 



£?sym 



' a b 

c c 

c c 

y b o d 



(4.38) 
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with the matrix elements given by, 



3 + cos N - 2 (2 X t) -2cos N - 1 {xt) 
8 

1 - cos jV - 2 (2 X t) 
8 ' 
N 2 - 4M 2 
4N(N - 1) ' 

3 + cos N ~ 2 (2 X t) + 2cos N ~ 1 (xt) 



4.2.2 Local invariants 

The two qubit local invariants Eq. (12.271 ) associated with the N qubit Kitagawa-Ueda state 
are given by 



T T 2,2,2 

J 3 = S S = Si + s 2 + s 3 

J 4 = s T Ts 



= ll 3 (l + co S (^ 2 )(2 X t)), 



-22- 3 co S 2 ^- 2 )( X t) sin 2 ( X t), 
^6 = £ijkSi{T s)j (T 2 s) k 



8 



2i = det T 




( ^ 2) (2xt)) 



2(7V " 2) (xt)) 



= 0. 



(4.39) 



We see that pairwise entanglement is manifest through the negative value of the invariant 
X5. Further, from Eq. (13.321) it is evident that X5 < implies that the state is spin squeezed. 




Figure 4.1: The invariants Z4, X4 — X| andXs, coiTesponding to a iV-qubit Kitagawa-Ueda 
state. Curve a : iV = 4, 6 : N = 6, and c : N = 8. 
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In Fig. 4.1, we have plotted the two qubit local invariants Z4, X4 — Zf and Z5, associated 
with the multiqubit kitagawa-Ueda state for different values of N. 



4.3 Atomic spin squeezed states 

Consider a system of N identical two-level atoms. The atomic squeezed states are denned 
by 

N \ 
J = -, Mj , (4.40) 



|^m) = A) exp (0 J 3 ) exp ( -i - J 2 



7T 



where 6 is conveniently parameterized as e 29 = tanh 2£. 

The atomic spin squeezed states \^m) are shown GUI to be the eigenvectors of the 
nonhermitian operator R3 , given by 



J_ cosh £ + J + sinh £ 
V / 2im¥2£ 



(4.41) 



where 



3 \*M 



Ml* 



at; 



(4.42) 



Here, we concentrate on the state 



|*o) = A) exp (9 J 3 ) exp ( -i — J 2 



7T 



N 

J = — , 
2 



(4.43) 



In the interaction of a collection of even number of iV two-level atoms with squeezed ra- 
diation, it has been shown [20] that the steady state |* ) of Eq. ( 14.431 ) is the pure atomic 
squeezed state for certain values of external field strength and detuning parameters. 
It is convenient to express the atomic squeezed state | \J/ ) as > 



Ac 



M=-J 



J ( K \ M8 
MO I — I e 



{AAA) 
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where the coefficients c^ (§ ) are given by l63l 



J\sJ{J + M)\ (J - M)\ 
2J 



J-M 




(J - M)!p! (p - M)\ (J + M - p)\ ' 



p=M 



(4.45) 



We now proceed to study the pairwise entanglement properties of atomic squeezed 
state |^o)- 

First and second order moments of the collective spin observable for the Atomic squeezed 
states: 

The first order expectation values of (Jj) are evaluated as follows: 

Consider the expectation value of the nonhermitian operator R3 (see Eq. (14.411 )) 



(tfol-Raltfo) 



1 



($o\J- cosh£ + J + sinh^ |*o) • 



V2 sinh 2£ 



(4.46) 



However, it is clear from Eq. (14.421) that 



(^ol^al^o) =0. 



So, we obtain, 



{^o\J- cosh£ + J+ sinh £ |^ ) = 0. 



(4.47) 



Further, it is evident that, 



<*ol4l*o> =0 



which in turn leads to, 



(* |^+eosh£ + J_ sinh £ |^ ) = 0. 



(4.48) 
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(J+) = 0, (J_)=0, 
or (Jx) = 0, (J 2 ) = 0. 

The first order moment ( J3) can be explicitly evaluated as follows: We have, 



J 3 |* > = io E Md J M0 {^)e^\JM), 

and therefore 

(^o|^3|^o) = 



M=-J 
J 



ai y: m 

M=-J 



MO 



TT 



(4.49) 



The second order moments ((JiJj + JjJi)) of the collective spin operator, are conve- 



niently evaluated in terms of the expectation values of the non-hermitian operator R 



2 . 



2 (J_ cosh £ + J + sinh £) ( J_ cosh £ + J + sinh £ ) 

3 ~~ 2 sinh 2£ 

J 2 cosh 2 £ + J| sinh 2 £ + (J_ J+ + J+J-) sinh £ cos £) 
~~ 2 sinh 2£ 

(^olfisl^o) = (J?- J 2 2 )coth2£ + (J 2 + J 2 )-i([J 1 ,J 2 ] + )coth2£. (4.50) 



since 

(^o|-R||*o) = 0, (see Eq. (14421 ) we obtain, 



(J 2 - J 2 2 ) coth 2£ + (J 2 + Jf ) - iflJi, J 2 ]+) coth 2£ = 0. (4.51) 



In other words, we have, 



Re ((#§}) = ( Jj 2 - J 2 2 ) coth 2£ + (J 2 + Jf } = 



(4.52) 



and 



Im (( J R 2 )) = (([J 1 ,J 2 ]+) coth 2£) = 0. 



'since J+ = and J_ = (J+) = 0, (J_) = (Ji) = (J 2 ) = 0. 



(4.53) 
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(J+ cosh£ + J_ sinh£)(J„ cosh£ + J + sinh£) 
2 sinh 2£ 

J + J_ cosh 2 £ + J_ J + sinh 2 £ + (J 2 + J 2 ) sinh £ cos £ 
2 sinh 2£ 

(Jf + J|) coth 2£ + ( Jl - J 2 2 ) + ( J 3 ) . (4.54) 

From Eq. ( 14.421 ), it is clear that (R^Rs) = 0, and therefore we get, 

(Jf + Jf> coth 2£ + ( J? - Jf) + ( J 3 ) = 
i.e., - (Jf - Jf) - (J 2 + J 2 ) coth 2£ = (J 3 ). (4.55) 

Simplifying the Eqs. (14.521 ). (14.55I ). we obtain the expectation values of J 2 and J 2 as, 

(Jl) = -\(Jz)e-* 

(jf> = -\(Js)e^. (4.56) 
Now, to compute the average value of J| , we use 

(J 2 ) = (Jf + Jl + Jl) = J(J + 1) (4-57) 

and obtain, 

<J 2 > = (J 2 -Jt-Jl) 

= J (J + 1) - ( J 3 ) cosh 2f . (4.58) 



-R3-R3 — 
<*0 1 4^3 1*0) = 



The second order expectation values (^o| [J\, ^l + l^o) an d (^ol [J2, Jsl + l^o) are deter- 
mined as follows: 

Let us consider ([J+, «/s]+) = (^o|<^+>^3|V'o) + (V'o|^3^+|'0o)- By computing each 
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term separately, 

J 

<* |J + J 3 |*o> = A Md J M0 (|) d J M , (|) e ( M+M ')V(^ + J - M' + 1) 5 M ,, M+ i 



M=-J 
J 



= 4, £ Md J MQ (|) <4 + io (|) eW+Wy/(J + M + 1)( J - AO 

Af=-J 

(J+J 3 ) = 0, (4.59) 



.a. 

Similarly, 



J 

'7T\ ,7 /VT 



M=-J 
J 



J 3 J+|*o> = 4> £ M4 (-) ( g ) e {M+M ' )e V(J ~ M')(J + WTT) 5 m >,m-i 



= 4> £ MdU (f ) <-io (f) eW+Dtyftj — M + 1)(J + M) 

M=-J 

(J 3 J+) = 0. (4.60) 
From Eqs. ( 14.591 (14.601 ). we obtain 

<[J+,J 3 ] + )=0. (4.61) 

So, we obtain, 

<[Ji,J 3 ]+> = 0, <[J 2 ,J 3 ]+) = 0. (4.62) 

We next determine the two qubit state parameters associated with the atomic squeezed sys- 
tems. 



2 The coefficient d J M0 (§) = V( J+M >j< J M >' , +A , 1 , A , (-1)^ for J + M = even and 

I 2 M 2 ^ - 

^mo(§) = f° r J -\- M = odd. Therefore, we obviously have d'l I0 d'[ I+1Q = 0. 
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4.3.1 Two qubit state parameters 

The components of the single qubit orientation vector s drawn from a collective atomic 
system, are given by (see Eqs. ( 14.491 ), ( 14.491 )) 



Sl = iv (Jl) = 0, 



S2 = ^2> = 0, 



S3 



M=-J 



%0 



7T 



3 (2M0) 



(4.63) 



Thus the average spin vector s for atomic squeezed states of Eq. ( 14.431 ) assumes the form 



0, 0, 



N 



The elements of the two qubit conelation matrix which are expressed in terms of the second 
order moments (see Eq. (13.111 )) are given by, 



t 



i.i 



1 

N - 1 



2 {^JiJj ~\- JjJij 

N 
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The diagonal elements of the correlation matrix T are obtained using Eqs. (14.561 ), (14.581 ) and 
are given by, 



tn 



*22 



*33 



1 



4[(Ji 2 )] 
N(N-l) N-l 



1 



N(N - 1) 



1 



(Js)e^-N 



-2(J 3 }e-^-N 
N(N-l) 

4[(J|)] 



1 



N(N - 1) N-l 



1 



N(N - 1) 
-2(J 3 )e 2 Z-N 
N(N - 1) 



N - 1 



N{N - 1) 

[J(J + 1) - ( J 3 ) cosh 2£ - iV] 

4(J 3 ) cosh(2Q + N 2 + JV 
N(N-1) 



(4.64) 



Further, from Eqs. (14.62b . (I4.53I ). it is easy to see that the off-diagonal elements of T are all 
zero 



tl2 


= *21 


= o, 


tl3 


= tai 


= 0, 


t 2 3 


= *32 


= 0. 



(4.65) 
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Thus, the correlation matrix T has the following structure, 



T = diag (ti,t 2 ,t 3 ) 



( -2(.h)e-^-N 

n(n-i) 






-2{J 3 )e 2 t-N 
N(N-l) 








4(J 3 ) cosh(2Q+jV 2 +jV 
N(N-l) 



\ 



4.3.2 Local invariants 



The two qubit local invariants (see Eq. (12.271 )) associated with the atomic spin squezed states 
are listed below: 



X 2 
X 3 

X 4 

X 5 



ti + ti + t 



3< 



T 2,2,2 
S S = S 1 + S 2 + s 3 

4 (J 3 ) 2 
N 2 ' 

s T Ts = slh + s 2 2 t 2 + a§ £ 3 



X 3 



4(J 3 ) cosh(2|£|) + iV 2 + iV 



N(N-1 

^ijk tlmn &i $1 tj m tk n 
2X 3 



N 2 (N 
e ijk Si (Ts)j (T 2 s) k 



w (2(J 3 )e~W+N) (2(J 3 )e^+N), 



0. 



(4.66) 



In Fig. 4.2, we have plotted the invariants X4, T4 — Xf and X5, as a function of the parameter 
x = e 2e ,for different values of N. These plots demonstrate that the invariant Z5 is negative, 
highlighting the pairwise entanglement (spin squeezing) of the atomic state. 




x 5 




-0.14 1 1 ' 1 ' 1 

0.2 0.4 0.6 0.8 1 

X 

Figure 4.2: The invariants I4, Z4 — Zf and2s, associated with the atomic squeezed state of 
AT two level atoms interacting with squeezed radiation. Curve a: N = A, b : N = 6, c : 

N = 8, and d : iV = 20. 
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In this chapter, we have considered few interesting symmetric multi-qubit dynamical models 
like Dicke states, Kitagawa-Ueda state generated by one axis twisting Hamiltonian and 
Atomic squeezed state |^o)- The density matrix of these states have a specific structure and 
belong to the special class of symmetric states Eq. ( 12.341 ) discussed in Chapter. |2] We have 
evaluated the two qubit local invariants and investigated the nonlocal properties associated 
with the states. In each case, the entanglement properties are reflected through the negative 
values of some of the two qubit invariants thus highlighting our separability criteria. 



Chapter 5 

Constraints on the variance matrix of 
entangled symmetric qubits 

The nonseparability constraints on the two qubit local invariants X4 < 0, X5 < and 
T4 — X| < derived, in Chapters |2] and serve only as sufficient condition for pair- 
wise entanglement in symmetric qubits. In this Chapter, we derive necessary and sufficient 
condition for entanglement in symmetric two qubit states by establishing an equivalence be- 
tween the Peres-Horodecki criterion B3l l44l and the negativity of the two qubit covariance 
matrix. Pairwise entangled symmetric multiqubit states necessarily obey these constraints. 
We also bring out a local invariant structure exhibited by these constraints. 

5.1 Peres-Horodecki inseparability criterion for CV states 

Peres-Horodecki inseparability criterion ll43l l44l viz., positivity under partial transpose 
(PPT) has been extremely fruitful in characterizing entanglement for finite dimensional sys- 
tems. It provides necessary and sufficient conditions for 2 x 2 and 2x3 dimensional systems. 
It is found that the PPT criterion is significant in the case of infinite dimensional bipartite 
Continuous Variable (CV) states too. An important advance came about through an identi- 
fication of how Peres-Horodecki criterion gets translated elegantly into the properties of the 
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second moments (uncertainties) of CV states Il42l . This results in restrictions Il42l |64l on 
the covariance matrix of an entangled bipartite CV state. In the special case of two-mode 
Gaussian states, where the basic entanglement properties are imbibed in the structure of its 
covariance matrix, the restrictions on the covariance matrix are found to be necessary and 
sufficient for inseparability Il42ll64ll . 

Here, we construct a two qubit variance matrix (analogous to that of CV states) for 
two qubits and derive corresponding inseparability constraints imposed on it. 

Let us first recapitulate succinctly the approach employed by Simon ll42l for bipartite 
CV states: The basic variables of bi-partite CV states are the conjugate quadratures of two 
field modes, 

i= (9i,Pl,92,jP2), (5-1) 
which satisfy the canonical commutation relations 

[ia,ip]=i^a/3, a, = 1,2,3,4 (5.2) 

where, 

n 

and J 

The matrix Q, is known as a symplectic matrix. All real linear transformations applied to 
the operators £ Q/ 3 that obey the commutation realtions Eq. (15.21 ), form a group known as the 
symplectic group. 

The second moments are embodied in the real symmetric 4x4 covariance matrix of 



J 
J 



1 
-1 



(5.3) 
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a bipartite CV state, which is defined through it's elements: 

^ = i({A| a ,A|g}), 
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(5.4) 



where, 



and 



(5.5) 



A| a , A^} = + A&A| a . (5.6) 

Under canonical transformations, the variables of the two-mode system transform as 

where 5 G Sp (4, R) corresponds to a real symplectic 4x4 matrix. Under such transfor- 
mations, the covariance matrix goes as 

V^V' = SVS T . 



It is convenient to caste the covariance matrix in a 2 x 2 block form: 



V 



( A 



C T B 



(5-7) 



The entanglement properties hidden in the covariance matrix V remain unaltered under a 
local Sp (2, R) (g) Sp (2, R) transformation. Such a local operation transforms the blocks 
A, B, C of the variance matrix Eq. (15.71 ) as 



A - 


-> A' 


= StASj, 




B - 


+ B' 


= S2 B S2, 




C - 


*C 


= SiCSf. 


(5.8) 
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There are four local invariants associated with V given in terms of the blocks A, B, C: 



det A 



det 5 



det C, 



Tr(AJCBC T J). 



(5.9) 



The Peres-Horodecki criterion imposes the restriction Il42l 



) 



2 



ii> 



(5.10) 



on the second moments of every separable CV state. 

The signature of the invariant Jg = det C has an important consequence: Gaussian states 
with /g > are necessarily separable, where as those with I' 3 < and violating Eq. (15.101) 
are entangled. 

In other words, for Gaussian states violation of the condition Eq. (15.10b is both necessary 
and sufficient for entanglement. 

In the next section, we explain a similar formalism for symmetric two qubits by con- 
structing a covariance matrix and analyzing its inseparability behaviour. 

5.2 Two qubit covariance matrix 

The basic variables of a two qubit system are expressed as a operator column(row) as 



The 6x6 real symmetric covariance matrix V of a two qubit system may be denned through, 




Oli,cr 2 j), i,j 



= 1,2,3. 



Vai;f3j — ^ ({A(ai,A(pj}) 



(5.11) 
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with a, (3 = 1, 2 ; i, j = 1, 2, 3. The variance matrix V can be conveniently written in the 
3x3 block form as 



/ 



V 



A C 
C T B 



(5.12) 



where, 



B 



2 [({^i,^}} - (o- 2j )] 

<5jj - (f2t) (0"2j) 



Hj r i r ji 



C 



- [(<ru(rzj) - (a u ) (a 2j )] 



(5.13) 



In other words, we have, 



A 
B 
C 



I — s s T , 



1 — rr 



T — sr 



(5.14) 



Here X denotes a 3 x 3 identity matrix and Sj, r{ and ty are the state parameters of an 
arbitrary two qubit state (see Eqs. (12.61 ) - (12.10b ). 

In the case of symmetric states, considerable simplicity ensues as a result of Eqs. ( 12.181 ), ( 12.191 ) 
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and the covariance matrix of Eq. ( 15.121 ) assumes the form: 



V 



A C 
C T A 



(5.15) 



where, 



Explicitly, 



and 



A = B = 1 - ss J 



and C = T - ss 1 . 



( I -si 



A 



SlS2 -S\S3 

SlS 2 1 - s\ -s 2 s 3 
\ -sis 3 s 2 s 3 I - si I 



(5.16) 



C 



fit - si 



tl2 - S1S2 t 2 2 - S 2 t 23 ~ S2S3 
tl3 ~ SlS 3 t 23 ~ S 2 S 3 t 33 - si 



(5.17) 



We now establish an important property exhibited by the off-diagonal block C of the 
covariance matrix of a symmetric two qubit state. 

5.3 Inseparability constraint on the covariance matrix 

Lemma: For every separable symmetric state, C = T — ss T is a positive semidefinite 
matrix. 
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Proof: Consider a separable symmetric state of two qubits 



P(sym-sep) = ^ Pw Pw ® Pw , ^Pw = 3-5 < p w < 1. 



(5.18) 



LU 



The state variables Sj and tij associated with a separable symmetric state have the following 
structure: 



Let us now evaluate the quadratic form n T (T — ss T ) n where n (n T ) denotes any arbitrary 
real three componental column (row), in a separable symmetric state: 



which has the structure (A 2 ) — (A) 2 and is therefore a positive semi-definite quantity.D 

This lemma establishes the fact that the off diagonal block C of the covariance matrix 
is necessarily positive semidefinite for separable symmetric states. And therefore, C < 
serves as a sufficient condition for inseparability in two-qubit symmetric states. 

We now investigate the inseparability constriant T — ss T < in the case of a pure 
entangled two qubit state. 




w 




(5.19) 



w 




(5.20) 
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An arbitrary pure two qubit state can be written in a Schmidt decomposed form 



|$) = Kl |0i0 2 ) + K 2 |li 1 2 ), kI + kI 



+ 4 = 1 



(5.21) 



where, 



< K2 < Ki < 1, 



are the Schmidt coefficients. The two qubit state can be written in the 4x4 matrix form 
using the basis {|0i 2 ), |0i 1 2 ), \U 2 ), |li 1 2 )} : 



I k\ 2kik 2 » 













(5.22) 



l 2kik 2 n\ J 



The 3x3 real symmetric correlation matrix T may be readily obtained using Eq. (12.101 ) 



as. 









-2«i« 2 , 
1 



TrT = 1. 



(5.23) 



The average qubit orientation (Eq. ( 12.91 )) has the form 



(0, 0, k\ - 4) = r. 



(5.24) 



From Eq. (15.231 ) and Eq. (15.241 ). it is clear that an arbitrary two qubit pure state is symmetric 
in the Schmidt basis. 
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The 3x3 matrix C = T — ss T takes the form 



C = T-ss i 










-2k\K2, 

A ,,1 ,.2 



4k|k| j 



(5.25) 



It can be clearly seen that C < 0, for all entangled pure two-qubit states. 

In other words, the condition C < is both necessary and sufficient for pure entangled 

two-qubit states. 

Interestingly, non-positivity of C completely characterizes inseparability in an arbi- 
trary symmetric two qubit state, which will be proved in the Sec. (I5.4I ). 



5.4 Complete characterization of inseparability in mixed two 
qubit symmetric states 



We prove the following theorem: 

Theorem: The off-diagonal block C of co variance matrix of an entangled two qubit mixed 
state is necessarily non-positive. 

Proof: An arbitrary two qubit symmetric state, characterized by the density ma- 
trix Eq. (12.171) . with the state parameters obeying the permutation symmetry requirements 
Eqs. ( 12.181 ), (12.191 ) has the following matrix form: 



Psym — ^ 



( 1 + 2 S 3 + *33 

A 
A 

\ (tll-t22) + 2tti2 



A* A* (tii-t 22 )-2iti2 

(t n +t 22 ) (in+t22) B* 
ihi + t 22 ) {t u +t 22 ) B* 



\ 



B 



B l-2s 3 + t 3 3 J 

(5.26) 
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in the standard two-qubit basis {|0i O2) , |0i I2) , 1 1 1 O2) , |li I2)} • 
Here, we have denoted, 

A = (si + i s 2 ) + (tia + **23) 
B = (si + is 2 ) - (ti3 + it 2 3)- 

Note that the two qubit basis is related to the total angular momentum basis | J, M) with 
J = 1, ; - J < M < J as follows: 



|0i0 2 ) 
|li 1 2 ) 
|0il 2 ) 

|li 2 ) 

and the following unitary matrix, 



IM>, 

|i, — 1 
1 



V2 



(|1,0) + |0,0)), 



-^(|1,0)-|0,0)), 



(5.27) 



(1 



u 







73 73 







V 



V2 



1 



73 °J 



(5.28) 



transforms the two qubit density matrix p sym of Eq. ( 15.261 ) to the angular momentum basis 
Eq. (15371 : 



U Psym U ] 



( 



PS 




where, 



Ps = - A 



( l + 2s 3 + i33 V2A* (t u -t 22 )-2it 12 ^ 

V2A 2(tu+t 22 ) V2B* 

y(t u -t 22 ) + 2it 12 V2B l-2s 3 +t 3 3 j 



(5.29) 
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So, an arbitrary two qubit symmetric state always gets restricted to the 3 dimensional 
maximal angular momentum subspace spanned by {| J max = 1, M) — 1 < M < 1} . How- 
ever, the partial transpose of p sym , does not get restricted to the symmetric subspace with 
J m£H = 1, when transformed to the total angular momentum basis Eq. (I5.27I ). 

Under the Partial transpose (PT) operation (say, on the second qubit), the Pauli spin 
matrices of second qubit change as 

C21 — > 021, 022 — > — (722, 023 — ► O23. 

When this PT operation, is followed by a local rotation about the 2-axis by an angle tt, the 
spin operators of the second qubit completely reverse their signs: 

Thus, PT map on the symmetric density operator Eq. (12.171 ) 



1 



3 3 



Psym = ^\ I ® I + Z^ s i + °2i) + CTli^jUj \ , (5.30) 

i=l i,j=l 



leads to 



1 ( 3 3 \ 

pjym = ^ I I ® -/" + ^{Vli Si ~ <T 2i Si) - CTij £T 2 j *ij j • (5-31) 
V i=l i=l / 

(Here T2 corresponds to partial transpose map on the second qubit). 
We thus obtain, 

( (hi + t 22 ) a* b* -(t n -t22) + 2it 12 \ 

a 1 + 2 s 3 + t 33 t 33 -l -a* 

b t 33 -l l-2s 3 + t 33 -b 

y -(til- 1 22 )- 2it 12 -a -b* (t n + t 22 ) J 

(5.32) ' 



t 2 .1 

Psym a 
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in the basis {|0i O2) , |0i I2) , 1 1 1 O2) , |li I2}} • Here we have denoted, 



a = -(si + is 2 ) - (ii3 + H23), 
b = (si+is2)-(ti3-it 2 3). 



Now a unitary transformation Eq. (15.271) which corresponds to a basis change Eq. (15.281 ) 

gives; 



pZ 2 m = u P T2 If 



{hi + t 22 ) 

->/2ti3 
-{hi - t 22 ) - 2iti2 



where, 



-V2ti3 -{tn - til) + 2iti2 a 

2t 33 b* 2 s 3 

b (in + 1 22 ) V2si 

2 s 3 V2si 4 



(5.33) 



a = V2{-si + is 2 + 2*23) 
6 = V2(*13 + «S2 + i*23)- 

It may therefore be seen that 3 dimensional subspace spanned by J max = 1 of total angular 
momentum of a symmetric two qubit state does not restrict itself to a 3 x 3 block form. 

Interestingly, a further change of basis defined by, 



\X) 



1 



7T 



|i,i>-|i,-i», 



\Y) = — (|1,1) + |1,-1)), 

\Z) = |i,o), 



(5.34) 
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which corresponds to a unitary transformation 



U' 



V2 



-10 10 

-i i 

y/2 

^ V2 J 



(5.35) 



on /)^y 2 m , leads to the following elegant structure 



-T' 2 
Psym 



'sym 



' ill *12 *13 Si 




Now a congruence ^| operation L pjy^ with 



(5.36) 



L 



X —s 



1 



gives, 



T n T ' 2 rt 

^ Psym 



in 



*12 
tl3 





*12 



*22 



^ 23 





iia 

*23 



^33 



si 



1 



1 Note that the congruence operation does not alter the positivity (negativity) of the eigenvalue structure of 
the matrix. 



Constraints on the variance matrix of entangled symmetric qubits 



87 



or 




1 / T-ss T 



2 



(5.37) 



V 







1 



It is therefore evident that negativity ofT — ss T necessarily implies negativity of the partially 
transposed arbitrary two qubit symmetric density matrix as, 



i.e., non-positivity of C = T — ss necessarily implies that partially transposed two qubit 
symmetric density matrix p^y m is negative. In other words, C < captures Peres 's insepa- 
rability criterion on symmetric two qubit state completely. Hence the theorem. □ 

In the next section, we explore how negativity of the matrix C reflects itself on the 
structure of the local invariants associated with the two qubit state. 

5.5 Local invariant structure 

The off-diagonal block C of the covariance matrix is a real 3x3 symmetric matrix and so, 
can be diagonalized by an orthogonal matrix O i.e., 



The orthogonal transformation corresponds to identical unitary transformation U ® U on 
the qubits. 

We denote the eigenvalues of the off-diagonal block C of the covariance matrix Eq. (15.121 ) 
by ci, C2, and C3. Restricting ourselves to identical local unitary transformations, we define 



pjym < ^ pj y 2 m < & Lf&tf < & C = T- ss T < 0; 



OCO^ = C d = (c u c 2 ,c 3 ). 
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three local invariants, which completely determine the eigenvalues c\, c 2 , C3 of C = T—ss T : 



Ti = det (C) = ci c 2 c 3 , 
J 2 = Tr (C) = ci + c 2 + c 3 , 

J 3 = Tr(C 2 ) = c 2 + c 2 + c 2. (5.38) 
The invariant Z 2 may be rewritten as 

T 2 = Tr(T - ss T ) = I -si, (5.39) 
since Tr (T) = 1 for a symmetric state. Here, we have denoted 

rp / T\ 2,2,2 2 

lr [s s ) = s 1 + s 2 + S3 = s . 

Another useful invariant, which is a combination of the invariants defined through 
Eq. (15.381 ), may be constructed as 

2 2 —2 

T 4 = = ci c 2 + c 2 c 3 + ci c 3 . (5.40) 



Positivity of the single qubit reduced density operator demands Sq < 1 and leads in turn to 
the observation that the invariant T 2 is positive for all symmetric states. Thus, all the three 
eigen values c\ , c 2 , C3 of C can never assume negative values for symmetric qubits and at 
most two of them can be negative. 

We consider three distinct cases encompassing all pairwise entangled symmetric states. 

Case (i): Let one of the eigenvalues c\ = and of the remaining two, let c 2 < and 
c 3 >0. 

Clearly, the invariant 2\ = in this case. But we have 



^4 = c 2 c 3 < 0, 



(5.41) 
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which leads to a local invariant condition for two-qubit entanglement. 

Case (ii): Suppose any two eigenvalues say, ci,C2, are negative and the third one C3 is 
positive. 

Obviously, T\ > in this case. But the invariant Z4 assumes negative value: 

24 = ci T 2 - c{ + c 2 c 3 < (5.42) 



as each term in the right hand side is negative. In other words, T4 < gives the criterion 
for bipartite entanglement in this case too. 

Case (iii): Let c\ < 0; c 2 and C3 be positive. 

In this case we have 

2i < 0, (5.43) 

giving the inseparability criterion in terms of a local invariant. 

The new set of local invariants (see Eqs. ( 15.381 ), ( 15.401 )) associated with the off-diagonal 
block C of the covariance matrix can be related to the symmetric two qubit local invariants 
given by Eq. (12.27b - In the following discussion, we restrict ourselves to the identical lo- 
cal unitary transformations U (8> U (Eq. ( 13.231 )) which transform the state vectors to the 
following form: 



(0,0, s ), 



and T 



( 



V «13 



1 23 



<-13 
t { ~ ] t" 



(5.44) 



t' 



33 / 
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det (C) = det (T - ss T ) 

= t~\_t ± (t' 33 - si) - (t23) 2 t~[ — (t'( S ) 2 t ± . 

= detT-s^jtJ. 

= ^i-f: (5-45) 



where we have used Eqs. (12.30b . (13.241) . We thus have, 



Ti=Jl-y. (5.46) 



We further find that, 

J 2 = Tr (C) = Tr (T — ss T ) 

= Tr (T) - Tr (ss T ) 

= 1 - s T s 

= 1-X 3 . (5.47) 

The invariant X3 , can be written as 

T 3 = Tr(C 2 ) = Ti[{T-ss T ) 2 } 

= Tr [(T 2 ) + (ss T ss T ) - Tss T - ss T T] 
= Tr (T 2 ) + (s T s) 2 - 2 s T Ts. 



From Eq. (12.27b . we have T3 given by, 



l 2 +l!-21 4 . (5.48) 



Now, we proceed to explore how this basic structure C < reflects itself via collective 
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5.6 Implications of C < in symmetric N qubit systems 

Collective observables are expressible in terms of total angular momentum operator as 



A' 



j=Y,\ & « (5 - 49) 



2 

where ff a denote the Pauli spin operator of the a th qubit. 

The collective correlation matrix involving first and second moments of J may be 
defined as, 

Using Eq. (13.71 ) and Eq. (13.81 ). we can express the first and second order moments (Jj), 
(JiJj + JjJi) in terms of the two qubits state parameters. After simplification, we obtain, 

( l_ tu (JV-l)t 12 (JV-1)* 13 \ 

{N-l)t 12 1 - t 22 + N (t 22 - s 2 2 ) {N-l)t 2 z 

\ (N - 1) t 13 (N - 1) t 23 l-t 33 + N (i 33 - a§) y 

(5.51) 



4 



We can now express in the following compact form, 



V W = ^ (J - SS T + (iV - 1) (T - SS T )) (5.52) 



where X is a 3 x 3 identity matrix; s T = [s\ , s 2 , s 3 ) and T denotes the two qubit correlation 

matrix (see Eq. (12.11I )). We simplify Eq. (15.521 ) further. 

By shifting the second term i.e., ^ ss T to the left hand side, we obtain, 

F W + ^ s tJ (i+(jv _i)c), C = T-ss t (5.53) 
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Expressing Sj in terms of the collective observables, we have, 

Si = — (JA = — Si. (5.54) 

N N 

Now, substituting Eq. (15.541 ) in Eq. (15.53b . we obtain, 

V W + }_SS T = j(l+(N -l)C). (5.55) 

For all symmetric separable states we have established that C > (see our theorem in 
Sec. |5.4| ). We thus obtain the following constraint on the collective correlation matrix yW : 

V (N) + 1 SS T <^1. (5.56) 
N 4 

Pairwise entangled symmetric multiqubit states necessarily satisfy the above condition. 

Note that under identical local unitary transformations U (8> U ® ... (8> U on the 
qubits, the variance matrix and the average spin 5 transform as 

y{N)> = OVWO T , 
and S' = OS, (5.57) 

where O is a 3 x 3 real orthogonal rotation matrix corresponding to the local unitary trans- 
formation U on all the qubits. Thus, the 3x3 real symmetric matrix + i SS T can 
always be diagonalized by a suitable identical local unitary transformation on all the qubits. 
In other words, (15.561 ) is a local invariant condition and it essentially implies: 

The symmetric N qubit system is pairwise entangled iff the least eigen value of the real 
symmetric matrix V^ N ' + SS T is less than N j '4. 
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5.7 Equivalence between the generalized spin squeezing inequal- 
ities and negativity of C 

Let us consider the generalized spin squeezing inequalities of Ref. |[3T1 

4(AJ|) 4(J fe P 

"Iv - < 1 ~ "/vs - ' (5 - 58) 

where Jf. = J • k, with k denoting an arbitrary unit vector, and 

(AJ|) = (J 2 k ) - (J k ) 2 . 

Expressing ( J k ) and ( Jf ) in terms of the two qubit state variables Sj and Uj we have, 

1 N 3 N 

< J *> = 9 E E ^ = t (s "*- k) (5 - 59) 



iV 1 ^ ^ 

= "J + 2 EE E ^^Pi) k i k J 
i,j a=l /3>a=l 

N 1 X ^ X ^ 

~~ "J + 2 ^ ^ j j 

i,j a,(3>a 

- N i N{N ~ 1] Tt kk 

i.e., (J|) = j (1 + (iV - 1) A; T TA:) . (5.60) 
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Using Eqs. ( 15.591 ), ( 15.601 ) in the generalized spin squeezing inequalities given in Eq. ( 15.581 ), 
we obtain, 



l[(j|)- ( j fc ) 2 ] < 



4 

N 



£ (l + (N-l)k T Tk)-^(s-kf 



< 1 



4 (N 2 (s-k) 2 
N 2 4 



(l + (N-l)k T Tk)-N(s-k) 2 < l-(s-k) 2 

{N-l)k T Tk < (N-l)(s-k) 2 

k T {T-ss T )k < 0, 
or C < 0. 



(5.61) 



Thus we find that the generalized spin squeezing inequality is equivalent to the condition 

C = T - ss T < 0. 
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We have constructed a two qubit variance matrix and have shown here that the off-diagonal 
block of the variance matrix C of a separable symmetric two qubit state is a positive semidef- 
inite quantity. An equivalence between the Peres-Horodecki criterion and the negativity of 
the covariance matrix C is established, showing that the covariance matrix criterion is both 
necessary and sufficient for entanglement in symmetric two qubit states. Thus symmetric 
two-qubit states satisfying the condition C < are identified as inseparable. Further, the 
inseparability constraint C < is shown to be equivalent to the recently proposed ll3Tj 
generalized spin squeezing inequalities for pairwise entanglement in symmetric iV-qubit 
states. An elegant local invariant structure exhibited by these constraints on the two qubit 
covariance matrix has also been discussed. 



Chapter 6 



Summary 

Quantum correlated multiqubit states offer promising possibilities in low-noise 
spectroscopy lfT8ll . high precision interferometry |[T6l [171 l65l and in the implementation 
of quantum information protocols |[66l . Multiqubit states which are symmetric under in- 
terchange of particles (qubits) form an important class due to their experimental signif- 
icance |[T3l l62l l67l as well as the mathematical simplicity and elegance associated with 
them. 

Individual qubits (two-level atoms) in a multiqubit system are not accessible in the 
macroscopic ensemble and therefore only collective measurements are feasible. Any char- 
acterization of entanglement requiring individual control of qubits cannot be experimentally 
implemented. For example, spin squeezing |fl9l , i.e., reduction of quantum fluctuations in 
one of the spin components orthogonal to the mean spin direction below the fundamental 
noise limit N/4 is an important collective signature of entanglement in symmetric N qubit 
systems and is a consequence of two-qubit pairwise entanglement lfT3ll30ll32l . 

Spin squeezing is one of the important quantifying signatures of quantum correla- 
tions in multiqubit systems. Spin squeezed atomic states are produced routinely in several 
laboratories (6l 13 today. However, it is important to realize that spin squeezing does not 
capture quantum correlations completely and it serves only as a sufficient condition for 
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pairwise entanglement in symmetric N qubit states. Investigations on other collective sig- 
natures ||3T1 l32l l33l of entanglement gain their significance in this context. 

In this thesis, we have investigated the pairwise entanglement properties of symmetric 
multiqubits obeying permutation symmetry by employing two qubit local invariants. We 
have shown that a subset of 6 invaraints {X\ — Xq}, of a more general set of 18 invariants 
proposed by Makhlin P6l . completely characterizes pairwise entanglement of the collec- 
tive state. For a specific case of symmetric two-qubit system, which is realized in several 
physically interesting examples like, even and odd spin states l49l . Kitagawa - Ueda state 
generated by one-axis twisting Hamiltonian |fT9l , Atomic spin squeezed states |[20l etc, a 
subset of three independent invariants is sufficient to characterize the non-local properties 
completely. For symmetric separable states, we have proved that the entanglement invari- 
ants Xi, Z4, T5 andXi — 2| assume non-negative values. 

Based on negative values of the invariants X±, X4, X5 and T4 — X|, we have proposed a 
detailed classification scheme, for pairwise entanglement in symmetric multiqubit system, 
Our scheme also relates appropriate collective non-classical features, which can be identi- 
fied in each case of pairwise entanglement. Further, we have expressed collective features 
of entanglement, such as spin squeezing, in terms of these invariants. More specifically, we 
have shown that a symmetric multi-qubit system is spin squeezed iff one of the entangle- 
ment invariant is negative. Moreover, our invariant criteria are shown to be related to the 
family of generalized spin squeezing inequalities ||3T1 (two qubit entanglement) involving 
collective first and second order moments of total angular momentum operator. 

Further, we have established an equivalence between the Peres-Horodecki criterion 
and the negativity of the off diagonal block C of the two qubit covariance matrix thereby 
showing that our condition is both necessary and sufficient for entanglement in symmetric 
two qubit states P4l . Pairwise entangled symmetric multiqubit states necessarily obey these 
constraints. We have also brought out an elegant local invariant structure exhibited by our 
constraints. 



Appendix A 

Pure and mixed density operators 



Consider a quantum system characterized by a state \ip) in the Hilbert space H. Using 
complete orthonormal basis {|ii n )} satisfying 

• orthonormality: (u m \u n ) = 5 mn 

• completeness: Yl n \u n ){u n \ = /, 

where I is the unit operator. 

We can expand \ip) as follows: 





(A. 



n 



The expansion coefficients c n satisfy the normalization condition 



i c n \ 2 = i Kiv> i 2 = 1- 



n 



n 
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Let us evaluate the expectation value of an observable A in the state 

(A) = (V|^|V) 

= y^y^(^in m )(M m |^|n w )(K n |-^) 

m n 

= ^2^2(u n \p\u m )(u m \A\u n ) 

m n 

— ^ ^ ^ ^ Pnm-^mn 
m n 

or (A) = Tr(pA), (A.2) 



where, 



P=|^)(VI (A-3) 



is the density operator associated with the quantum system (elements of which are denoted 

by (u n \p\u m ) = p nm in the given basis {\u n )}). 

The density operator satisfies the following properties: 



1 . p is a hermitian matrix 



2. p is positive semi definite, 



3. p has unit trace, 



Pmn ~ Prim- (A.4) 



P > 0. (A.5) 



Trp = l. (A.6) 
From the structure of p = (see Eq. (1A.3I) ) it is clear that, 



p 2 = p, (A.7) 



i.e., p is idempotent operator. 
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Density operators satisfying Eq. dA.31 ) form a subclass of the more general class of 
density operators and are termed as pure density operators in contrast to mixed states. 

Mixed density operators are, in general, a convex mixture of pure states : 



(A.8) 



Number of real independent parameters characterizing the density operators: 

Writing a pure density operator p (IA.3b in the basis explicitly, we obtain, 



' | C! | 2 C\C* 2 

c 2 c* |c 2 | 2 



• ClC, 

• c 2 < 



V 



\Cn\ 2 ) 



(A.9) 



i,.e, the density matrix given above is specified completely by the complex coefficients c n , 
which are constrained by the normalization condition ^ n | c n | 2 = 1. In other words, the 
density matrix of Eq. (IA.3I ) is completely characterized by (2n — 1) real parameters. 

However, the more general class of mixed density operators satisfying the properties 
(Eqs. (IA.4b - (1A.6I) ) are characterized by n 2 — 1 real parameters as shown below: 
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The density operator p of a mixed quantum state may be explicitly written in the matrix 
form (in a suitable complete orthonormal basis {|« n )}) as > 



Pn P12 

Pl2 P22 



Pin 
P2n 



\ 



Tr (P) = ^2 Pnn = I- 



(A. 10) 



\ Pnl Pn2 ■ ■ ■ Pnn J 

We count the number of parameters as follows: 

• diagonal elements of p are real and are constrained by the unit trace condition. So, 
we have (n — 1) real parameters specifying the diagonal elements. 

• There are n ^ n ~ 1 ^ independent complex parameters p nm = p* mn ->n ^ in, which fix 
the off-diagonal elements of p. This leads to n(n — 1) real parameters. 

So, the total number of independent real parameters characterizing a mixed quantum 
state p are given by (n — 1) + n(n — 1) = n 2 — 1. 

Single qubit density matrices: 

There are three real parameters specifying a mixed density operator of single qubit system 
which is represented by a 2 x 2 hermitian matrix of unit trace as, 



(A. 11) 




The Pauli spin matrices, 



1 



1 



-i 



i 



1 

-1 



(A. 12) 
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together with the identity matrix 

I 



/ 



V 



1 
1 



(A. 13) 



provide a matrix basis for any 2x2 matrices. 

In the case of spin-i quantum states, we can always express 



p = -[I + erisi + <J 2 S2 + a s s 3 ], (A. 14) 



where, it is readily seen that 

si = (o"i) = Tr (pax), 
S2 = (o" 2 ) = Tr (pa 2 ), 

s 3 = (a 3 ) =Tr(pa 3 ). (A. 15) 

Two qubit quantum system: 

There are n 2 — 1 = 15 parameters characterizing a two qubit system. A natural opera- 
tor basis {/ (g) I, an, a 2 i, <Jn(T2j} is employed generally to expand the two qubit density 
matrix: 

1 I 3 3 3 

UjCFli02j I , (A. 16) 



j / 3 3 3 

p = - 1 1 (8> i + ^2 Si an + ^2 ^ n + U 

\ i=l i=l i,i=i 



with 



Oli = 

a 2i = I®Gi. (A. 17) 



The single qubit state parameters Sj and are given by 



si = Tr (pan) 

n = Tr (pa 2i ), i = 1,2,3, (A.18) 
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and the remaining two qubit correlation parameters are evaluated as follows: 

t tj =Tr(pa u a 2j ). (A. 19) 



Appendix B 

Peres PPT criterion 



Any bipartite state p defined on the Hilbert space H is separable if it can be expressed in 
the convex product form 



Here, p w and p w are the density operators defining the systems 1 and 2 respectively. Quan- 
tum systems which cannot be written in the form Eq. (IB. lb are called entangled states. The 
fundamental question in entanglement theory is the following: Given a composite quantum 
state how do we test if the state is entangled? 

In this context, Peres B3l has shown that the partially transposed density matrix of 
a separable state is a positive definite matrix. In other words, negative eigenvalues of a 
partially transposed density matrix necessarily imply entanglement in a quantum state. 

To demonstrate this explicitly, let us write the matrix elements of the separable state 
given in Eq. ( IB. 11 ) in the standard basis. 



Psep = y^^PwPh' ® PW where < p w < 1 and = 1. 



(B.l) 



w 



111 



(mp\p sep \nv) 



m 



p$\n) (n\p$W) 



w 



or 




(B.2) 



w 



104 



Peres PPT criterion 105 

Here, Latin indices refer to the system 1 and Greek indices correspond to the system 2. 

The partial transposition (PT) of any density matrix p with respect to one of the sub- 
systems, (say system 2) is defined as, 



Pmpu\nv ~ Pmvw 



Thus, the partial transpose of p scp (Eq. (1B.2I )) is given by 

= ^m^w^w- (B - 3) 

to 

Using the hermiticity property, p\ = pi (see Eq. (1A.4I )) of the density matrices, we obtain, 

(/i 2) ) T = (/i 2) r, (b.4) 

and thus, 

{Psep)m^\nu = ^ ] Pw (p w )mn (Pv) )i>n 
w 

or P^ p = 52pu,(p$)® (J®)*. (B.5) 

w 

Since correspond to physical density matrices, p^p i s again a physically valid 

separable density matrix. 

Therefore, the PT operation preserves the trace, hermiticity and also the positive semi 
defmiteness of a separable state whereas the last property need not be respected by an 
entangled state. 

Horodecki et. al P4l showed that Peres' positivity under partial transpose (PPT) is 
both necessary and sufficient for entanglement in 2 x 2 and 2x3 systems. However, 
for higher dimensions, there exist bound entangled states which are PPT states, though 
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they are inseparable. Negative under partial transpose (NPT) is a sufficient condition for 
entanglement in higher dimensional composite quantum systems. 



Appendix C 

A Complete set of 18 invariants for 
an arbitrary two qubit state 

An arbitrary two qubit density matrix (see Eq. (I2.6I )). specified by 15 real parameters 
{si, rj, tij}, is characterized by a complete set of 1 8 polynomial invariants ll36l . These local 
invariants are given in terms of the state parameters associated with an arbitrary two qubit 
system. Any two density matrices p\ and p2 are said to be locally equivalent if and only 
if all the 18 invariants (see Table. 1) have identical values for these states. To illustrate that 
these 18 invariants form a complete set, we may chose to work in a basis in which the 3x3 
real correlation matrix T (which is nonsymmetric, in general) is diagonal. Such a singular 
value decomposition of T can be achieved by proper rotations G 50(3, R) : 

T d = (1) ro (2)T = diag(tij ^ h) 
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We may note here that the diagonal elements of the correlation matrix T d are not the eigen- 
values of T as T 0^ T is not similarity transformation. However, 



and 



rpd rpcf^ 



Q(l) T q{2)T q{2) jiTQfl) 1 

(D rT T (ir 



t\ 
t\ 



t 2 



3 / 



rpd T rpd _ q{2) rpT rp 



(C.2) 



In other words, (t\, t?,, £§) are the eigenvalues of real symmetric matrix TT T as well as 
T T T. In order to determine the eigenvalues (t 2 , t 2 ., t 2 ) the following polynomial quantities 
may be employed: 



det (T T ) 
Tr (T T T) 
Tr [(T T T) 2 ] 



t\ + t\ + t\. 



(C.3) 



Note that det (T T T), Tr (T T T), Tr (T T T) 2 are invariant under local unitary transforma- 
tions on the two qubit state. It is easy to see that 



det (TT T ) = det (T)det (T T ) = (det (T)) 2 , 

= det(0Wr*0( 2 ) T )det(0( 2 ):K, o(D T ), 

= det(r d )det(T d ), 

= (det(T d )) 2 . 



(C.4) 
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Thus | det(T) | itself may be used in the first line of Eq. dC3b instead of det (TT T ). So, 
Makhlin chooses the first three elements of his set of local invariants as, 

h = det(T) =txt 2 t 3 , (C.5) 
I 2 = Ti{T T T)=t\ + tl + tl (C.6) 
h = Tr[(T T T) 2 ]=tf + 4 + 4 (C.7) 

The diagonal form of T viz., (tx, t 2 , £3) can be determined using the invariants Ix-3 up to 
a simultaneous sign change for any two of them. Now, a local rotation on the first qubit 
R(tt) 1 <8> I (where R(tt) 1 is the it rotation about the axis i = 1,2, 3), may be used to fix the 
signs of tx, t 2 , £3. It is then convenient to adopt the convention, (i) if Ix > 0, elements of 
T d are all positive and (ii) tx, t 2 , t%, are all negative, when Ix < 0. 

Let us restrict to two qubit states with a fixed diagonal correlation matrix T ( which is 
achieved through appropriate local unitary operations on the qubits). 

To determine the absolute values of the state parameters, sx, s 2 , S3, the invariants 
of Table 12. ll are used: 

T T 2,2,2 

14 = s s = s 1 + s 2 + S3, 

h = s T TT T s = s\ti + s 2 2 tl + sltl, 

I 6 = s T (TT T ) 2 s = s\t\ + s%4 + s% 4 (C- 8 ) 

(Here, T is considered to be nondegenerate, i.e., tx 7^ t 2 7^ t%). The absolute values of state 
parameters n, r 2 , rs, can be determined via the invariants I-j^g of Table |2~T1 

T T 2,2,2 

7 7 = r r = r\ + r 2 + r 3 , 

7 8 = r T TT T r = r 1 2 t2 +r 2 t 2 + r 2 t 2 ) 

J 9 = r T (rr T ) 2 r = r 1 2 ^ + r 2 2 t| + r 3 2 t|. (C.9) 
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Now, the invariant, Iiq (hi) is useful to fix the overall sign of s±, s 2 , and S3 (n, r 2 , and 

^10 = e ijk Si(TT T s)j ([TT T ] 2 s) k 

= {4 (4 - 4) + 4 (tj - 4) + 4 (4 - 4)) Sl S2 S3 , 

k 

4 (4 - 4) + 4 (4 - 4) + 4 (4 - 4)) n r 2 r 3 . (c.io) 



hi = e ljfe r i (T T Tr) J ([r T r] 2 r) 



Furthermore, the relative signs between s^r, are determined using the invariants /12-14, 



112 = s T Tr = sir 1 ti + s 2 r 2 t 2 + s 3 r 3 t 3 , 

1 13 = s T TT T Tr = 3^4 + 82^4 + ^4, 

hi = eijkQm.nSiritj m t kn = sirit 2 t 3 + s 2 r 2 tit 3 + s 3 r 3 t\t 2 , (C.ll) 



which provide three linear constraints on sin., s 2 r 2 and s 3 r 3 . 

Next, the individual signs of (si, s 2 , s 3 ) and (ri, r 2 , r 3 ) can also be determined 
when, atleast two components, say, si, s 2 are nonzero, and 53 = 0. In this case, the signs of 
si and s 2 can be made positive with the help of local rotations^] (i?(7r) 1 (g)i?(7r) 1 ), (R(tt) 2 (>S> 
R(it) 2 ). So, with si, s 2 > and S3 = the invariant /15 of Table [27T1 has the form 

h 5 = s lS2 t 3 r 3 [4-4], (C.12) 

thus fixing the sign of r 3 , provided £3 ^ 0. If t 3 = 0, the invariant In ( which is evaluated 

for si, s 2 > and S3 = £3 = 0) 

In = s 1 s 2 t 1 t 2 r 3 [t 2 2 - 4] (C.13) 



is utilized for determining the sign of r 3 . 

'Note that the diagonal form of T remains unchanged under R(i^Y ® R^T with i = 1, 2, 3. 
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A similar argument is applicable using - the invariants /11,16,1s f° r fixing the sign of 
S3 (^3 7^ 0)> when r\ and T2 are nonzero and = 0. 
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